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(57) Abstract 

A method and an apparatus measure blood oxygenation in a subject. A first signal source applies a first input signal during a first 
■; time interval. A second signal source applies a second input signal during a second time interval. A detector detects a first parametric 
signal responsive to the first input signal passing through a portion of the subject having blood therein. The detector also detects a second 
parametric signal responsive to the second input signal passing through the portion of the subject. The detector generates a detector output 
, signal responsive to the first and second parametric signals. A signal processor receives the detector output signal and demodulates the 
detector output signal by applying a first demodulation signal to a sign^ responsive to the detector ouQ)ut signal to generate a first output 
signal responsive to the first parametric signal. The signal processor applies a second demodulation signal to the signal responsive to the 
detector output signal to generate a second output signal lesponsiveto the second parametric signal. The first demodulation signal and the 
second demodulation signal both include at least a first component having a first frequency and a first amplitude and a second compcment 
having a second frequency and a second amplitude. The second frequency is a harmonic of the first frequency. The second amplitude is 
related to the first amplitude to minimize crosstalk from the first parametric signal to the second output signal and to minimize crosstalk 
from the second parametric signal to the first ou^ut signal. 
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Field of the Invention 

The present invention relates to the field of signal processing, and, more 

particularly, relates to the field of processing of signals generated in a physiological 
15 monitoring system, such as, for example, in a system for measuring blood oxygen 

saturation using pulse oximetry* 

Description of the Related Art 

The present invention will be described herein in connection with a pulse oximetry 

apparatus and a method, which are used to measure blood oxygen saturation in a subject, 
20 such as, for example, a human patient The teachings of the present invention can be used 

in other applications wherein useable signal information is obtained ia a noisy 

environment 

In an exemplary pulse oximetry apparatus and a corresponding method, blood 
oxygen saturation is determined by transmitting pulses of electromagnetic energy throug}i 

25 a portion of a subject which has blood flowing therein (e.g., through a finger, through an 
ear lobe, or other portion of the body where blood flows close to the skin). In the 
examples described herein, the pulses of electromagnetic energy comprise periodic pulses 
of red light having waveloigths of approximately 660 nanometers, for example, and 
periodic pulses of infirared light having wavelengths of approximately 90S nanometers. As 

30 described, for example, in US. Patent No. 5,482,036 and in U.S. Patent No. 5,490,505, the 
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pulses of red light and the pulses of infiared li^t are applied with the same poiodicity but 
in an alternating and non-overlapping manner. In particular, in preferred embodiments, the 
red pulses are active for approximately 25% of each cycle and the infrared pulses are also 
active for approximately 25% of each cycle. The red pulses are separated in time from the 
infrared pulses such that both pulses are inactive for approximately 25% of each cycle 
between a red pulse and the next infrared pulse and both pulses are inactive for 
approximately 25% of each cycle betwem an infrared pulse and the next red pulse. 
(Although described herein below in connection with pulses having 25% duty cycles, it 
should be understood by persons of skill in the art that the duty cycles of die pulses can be 
changed in some applications.) A&er propagatmg through the portion of the subject, the 
red pulses and the infrared pulses are detected by a detector ^ch is responsive to li^t at 
both wavelengths and which goierates an electrical signal whidi has a predictable 
relationship to the intensity of the electromagnetic energy incident on the detector. The 
electrical signal is processed in accordance with the present invention to provide a 
representation of the blood oxygen saturation of the subject 

In conventional time division multiplocing (TDM) demodulation that uses 
rectangular waves to drive the red and infrared LEDs, the conventional process of 
demodulation using square waves can result in the aliasing of the ambient noise 
components that come close to the sidebands of harmonics and the frmdamental frequency 
of the rectangular waves, and the noise components are thus collapsed into the output 
signal generated by the demodulation. In particular, it is very difficult to avoid including 
hatmonics of the line frequency in the demodulated output signal. 



The present inv^on avoids the problems associated with conventional 
demodulation and separation of TDM signals. In particular, the presot invention avoids 
the problem of aliasing of die ambient noise into the passband of the system by selectively 
demodulating certain harmonics of the TDM signal. For example, in one embodiment, 
only two harmonics (e.g., the frmdamental and the first hannonic) are demodulated. In 
other embodiments, more hannonics axe demodulated. The present invention specifically 
addresses solutions to problems caused by crosstalk resulting from filtering and also 



Summarv of the hvention 
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resulting fiom demodulating with only ceitain hannonics instead of demodulating with all 
harmonics as is done using conventional square wave demodulation. In a digital 
implementation of the present invention, the ou^ut of the photodetector is initially 
sampled at a very high frequency (e.g., 46,875 Hz), and the signals are decimated ( where 
decimation is lowpass filtering foUowed by sample rate compression) such that the final 
output signals arc generated at a relatively low sampling rate (e.g., 62.5 Hz) v^di 
provides increased resolution at the output. Thus, bandwidth is traded for resolution in the 
output signal, thus increasing the signal to noise ratio. 

One aspect of the present invention is an apparatus for measuring blood 
oxygenation in a subject The apparatus conptises a first signal source whidi applies a. 
first mput signal during a first time interval. A second signal source appUes a second mput 
signal during a second time interval. A detector detects a first paramedic signal respcmsive 
to the first input signal passing through a portion of ttie subject having blood thoein. The 
detector also detects a second parametric signal responsive to the second input signal 
passing through the portion of the subject The detector generates a detector output signal 
responsive to the first and second parametric signals. A signal processor receives flie 
detector output signal. The signal processor danodulates the detector outpvX signal by 
applying a first demodulation signal to a signal-responsive to the detector ou^ut signal to 
genraate a first output signal responsive to the first parametric signal and by applying a 
second donodulation signal to the signal req)onsive to the detector output signal to 
gaiciate a second ou^ut signal responsive to Ac second parametric signal. Each of the 
first demodulation signal and flie second demodulation signal conqnises at least a first 
component having a first fiwquency and a first amplitude and a second component having a 
second frequency and a second amplitude. The second fixquency is a harmonic of the first 
frequency. TIjc second ampKtiide is selected to be related to the first ampUtude to 
minimize crosstalk fiwm the first parametric signal to the second output signal and to 
minimize crosstalk from the second parametric signal to the first output signal. In one 
embodiment the second amplitude is determined by turning off one of the first and second 
signal sources and measuring the crosstalk between one of the parametric signals and flie 
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non-coiresponding output signal while varying the second amplitude. A second amplitude 
is selected that minimizes the measured crosstalk. 

Another aspect of the present invention is a method of minimizing crosstalk 
between two signals generated by applying a first pulse and a second pulse to measure a 
parameter. The first pulse and the second pulse are applied periodically at a first repetition 
rate defining a period The first pulse is generated during a first interval in each period, 
and the second pulse is genoated during a second interval in each period. The second 
interval is spaced apart fiom the first interval. The first and second pulses produce first 
and second parametric signals responsive to the parameto'. The first and second 
parametric signals are received by a single detector that ou^uts a composite signal 
responsive to the first and second parametric signals. The method comprises the stq> of 
applying a first demodulation signal to the composite signal to goierate a first 
demodulated ou^ut signal wherein the first draiodulation signal comprises at least a first 
component having a first fiequency corresponding to the first repetition rate. The first 
component has a first amplitude. The first demodulation signal fiirther comprises a second 
component having a second firequency that is a harmonic of the first firequency. The 
second component has a second anq)litude which has a selected proportional relationship 
to the first amplitude. The method fiirther includes the step of applying a second 
demodulation signal to the composite signal to generate a second demodulated output 
signal. The second demodulation signal conqmses the first con^onent at the first 
frequency and the first amplitude and further comprises the second componrat at the 
second frequency and the second amplitude. At least one of the first and second 
components of the second demodulation signal has a selected phase difference with respect 
to die corresppnding one of the first and second components of the first demoduladon 
signal. The method fiirther mcludes the steps of lowpass filtering the first demodulated 
output signal to generate a first recovered output signal responsive to the first parametric 
signal; and lowpass filtering the second demodidated output signal to generate a second 
recovered output signal responsive to the second parametric signal 



Substitute Sheet (Rule 26) 



o ■ o 

wo 99/52420 PCT/US99/07825 



Preferably, the selected phase difToieiice is Also preferably,- the first pulse and 
the second pulse are generally rectangular pulses having a respective duty cycle. The 
rectangular pulses comprise a plurality of sinusoidal components including a fundamental 
component corresponding to the first frequency and a first harmonic component 
5 corresponding to the second firequency. The fimdamental component has a fimdamental 
component amplitude and the first hamionic component has a first harmonic component 
amplitude. The first harmonic component amplitude is related to the fimdamental 
harmonic component amplitude by a first proportionality value. The second amplitude of 
the second component of the first demodulation signal is related to the first amplitude of 
10 the first component of the first dOTiodulation signal by a second proportionality value 
which is approximately the inverse of the first proportionality value. 

The method in accordance with this aspect of the invention preferably includes the 
fijrther steps of sampling the composite signal when neither the first pulse nor the second 
pulse is active to obtain a sampled signal; and measuring the sampled signal to detemiine a 
15 noise level of the parametric signals. 

In a fiirther embodiment according to this aspect of the present invention, the 
method fiirther includes the steps of performing a transform on the composite signal to 
generate a spectra of the composite signal; sampling the spectra at a plurality of 
fi^uencies other than at predetenoMned ranges of frequencies around the first firequency 
20 and around harmonics of the first frequency; determining an average of the magnitudes of 
the sampled plurality of fitsquencies; and comparing the average to a selected threshold to 
dctennine whether the average magnitude exceeds the selected threshold. 

Another aspect of the present invention is a method of demodulating a composite 
signal generated by applying first and second periodic pulses of electromagnetic energy to 
25 a systOTi having a parameter to be measured and by receiving signals responsive to the 
electromagnetic energy after having passed through the system and being affected by the 
parameter being measured. The signals are received as a composite signal having 
components responsive to the first and second pulses. The method comprises the stiep of 
qiplying a first demodulation signal to the composite signal to generate a first 
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demodulated signaL The first demodulation signal comprises a fitst component having a 
first firequency coirespbnding to a repetition firequency of the first and second pulses and 
comprises a second component having a fiequency that is a hannonic of the fiist 
fiequency. The first component has a first amplitude and the second component has a 

5 second amplitude. The second amplitude has a predetermined relationship to the first 
anq)iitude. The predetermined relationship is selected to cause the first demodulated 
signal to have low firequency components responsive only to the first pulse. The method 
includes the fitrther stq) of lowpass filtering the first demodulated signal to genmte a first 
output signal. The first output signal varies in response to an effect of the parameter on the 

10 electromagnetic energy received fi:om the first pulse. 

Preferably, the method in accordance with this aspect of the invention includes the 
fijrther step of applying a second donodulation signal to the composite signal to generate a 
second demodulated signal. The second demodulation signal has first and second 
components corresponding to the first and second con^nents of the first demodulation 

15 signal. At least one of the first and second components of the second demodulation signal 
has a selected phase relationship with the corresponding one of the first and second 
components of the first demodulation signal. The method includes the further step of 
lowpass filtering the second demodulated signal to generate a second output signal. The 
second output signal varies in response to an effect of the parameter on the electromagnetic 

20 energy received from the second pulse. 

Another aspect of the present invention is a pulse oxim^iy system that comprises a 
modulation signal generator. The modulation signal generator gmerates a first modulation 
signal that comprises a first pulse that repeats at a first repetition fipequenQT. Thefirstpulse 
has a duty cycle of less than 50%. The modulation signal gmeiator gmerates a second 

25 modulation signal comprising a second pulse that also repeats at the first repetition 
frequency. Thesecondpulsehasaduty cycle of less than 50%. The second pulse occurs 
at non-overlapping times witih respect to the first pulse. Each of the first and second pulses 
comprises a plurality of components wherein a first component has a firequency 
corresponding to the rq)etition fi^uency and wherein a second component has a second 

30 fiequency corresponding to twice the first firequency. The second component has an 
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amplitude which has a first predetennined relationship to an amplitude of the first 
component. A first transmitter emits electromagnetic energy at a first wavelength in 
response to the first pulse; and a second transmitter emits electromagnetic energy at a 
second wavelength in response to the second pulse. A detector receives electromagnetic 
energy at the first and second wavelengths after passing through a portion of a subject and 
generates a detector output signal responsive to the received electromagnetic energy. The 
detector output signal includes a signal component responsive to attenuation of the 
electromagnetic energy at the first wavelength and a signal component responsive to 
attenuation of the electromagnetic energy at the second wavelength. A first demodulator 
multiplies the detector signal by a first demodulation signal and generates a first 
demodulated output signal. The first donodulation signal comprises a first component 
having the first frequency and having a first amplitude. The first demodulation signal also 
comprises a second component having the second firequency and having a second 
amplitude. The second amplitude has a second predetermined relationship to the first 
amplitude. The second predetermined relationship is approximately inversely proportional 
to the first predetermined relationship. A second demodulator multiplies the detector 
signal by a second demodulation signal and generates a second demodulated output signal. 
The second demodulation signal comprises a first component having the first firequency 
and having tixe first amplitude. The second demodulation signal fiirther comprises a 
second component having the second fi*equency and having the second amplitude. At least 
one component of the second demodulation signal has a selected phase relationship with a 
corresponding one component of the first demodulation signal Preferably, the selected 
phase relationship is a n phase difference. 

Another anbodiment incorporates decimation before demodulation. In yet another 
embodimOTt, a multi-channel demodulator, with or without pre-demodulation decimation 
is disclosed. 

In yet another embodiment, an adaptive algorithm is used to control the operation 
of pre-demodulation decimators and post-demodulation decimators. The adaptive 

-7- 



Substitute Sheet (Rule 26) 



o 

wo 99/52420 



PCTAJS99/07825 



algorithm may control both the chazacteiistics of a lowpass filter in the decimator and the 
decimation rate provided by a signal rate compressor in the decimator. 

Another embodiment of the invention is a method for selecting a sample rate that 
reduces the interference caused by ambient lig^t 

Brief Description of the E>rawings 

The present invention will be described below in connection widi the 
accompanying drawing figures in which: 

Figure I illustrates an exemplary block diagram of a representation of a signal 
processing system in accordance with 'the present invention used to determine blood 
oxygen saturation in a subject; 

Figure 2 illustrates aemplary waveforms of the current throu^ the LEDs in 
Figure 1 and the resulting intensities of the red lighted the infirared light generated by the 
LEDs; 

Figure 3 illustrates a block diagram of the overall processing system in 
accordance with the present invention; 

Figure 4 illustrates a frequency spectra of the first modulation signal M\(t) for n = 
0, 1, 2, • . where the horizontal axis represents frequency and the vertical axis represents 
the energy in the DC and harmonic components of the signal; 

Figure 5 illustrates an exemplary spectrum of the first and second harmonics of the 
present invention when the fimdamental fiiequmcy is selected to be 316.7 Hz in 
comparison to the fimdamental and harmonics of conventional 60 Hz power; 

Figure 6 illustrates the effecc of the value of B on tiie measured signal ovxpvi 

Sj {t) responsive to the red modulation pulses as the value of B is varied while the infrared 
modulation pulses are off; 

Figure 7 illustrates a preferred embodiment of the present invention implemented 
in a digital processing system; 

Figure 8 illustrates a detailed block diagram of the demodulation portion of the 
present invention; 

-8- 



Substitute Sheet (Rule 26) 



wo 99/52420 

4 



PCT/US99/07825 



Figure 9 illustrates a detailed block diagram of the modulation portion of the 
present invention; 

Figure 10 illustrates the red drive waveform and the infrared drive waveform 
generated by the modulation portion of Figure 9; 

Figure 1 1 illustrates the demodulation waveforms generated by the demodulation 
portion of Figure 8; 

Figure 12 illustrates a method of time domain sampling the digital detection signal 
during the times when both the red pulses and the infrared pulses are off to obtain 
information regarding the level of ambient noise; 

Figure 13 illustrates a block diagram of a system that performs the time domain 
sampling of Figure 12; 

Figure 14 illustrates a method of frequency donudn sampling to determine the 
noise floor at frequencies other than the signal fi^uencies; 

Figure 15 illustrates a block diagram of a system that performs the frequency 
domain sampling of Figure 14; 

Figure 16 illustrates a block diagram of the overall processing system in 
accordance with a pre-demodulation decimation embodiment of the present invention; 

Figure 17 illustrates a block diagram of a multi-channel processing system in 
accordance with a pre-demodulation decimation embodiment of the present invention; 

Figure 18 illustrates a block diagram of an adaptive multi-channel processing 
systOTi in accordance with a pre-demodulation decimation embodiment of the present 
invention; 

Figure 19 illustrates a flowchart of a method for choosing the modulation 
frequency and decimation rate in order to minimize the affects of ambient light; and 

Figure 20 is a gr^h to be used in connection with graphical method for designing 
a demodulation system to minimize interference due to ambient ligjbt 
Detailed Description of the Preferred Embodiment 

Figure 1 illustrates an exemplary block diagram of a representation of a signal 
processing system 100 in accordance widi the present invention used to determine blood 
oxygen saturation in a subject, such as, for example, a human subject- In the example 
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presented, the measurements are performed on a portion of the subject, such as a finger 
102 illustrated in Figure L An LED modulation circuit 104 drives a pair of back-to-back 
light emitting diodes (LEDs) 106, 108 by applying a periodic signal to the two light 
emitting diodes 106, 108. The LED 106 is selected to emit electromagnetic energy in the 
red visible light range, and has a wavelength of, for example, approximately 660 
nanometers. The LED 108 is selected to emit electromagnetic energy in the infiared 
range, and has a wavelength of^ for example, approximately 905 nanometers. The T |yp 
modulation circuit 104 suppUes current in alternating directions so that the two LEDs 106, 
108 are activated one at a time. In particular, as illustrated by a current waveform 120 in 
Figure 2, current is first applied in a forward direction with respect to the red T.TTp io6 
during a first time intenral 122 having a duration x. Thereafter, no current is applied to 
either LED during a second time interval 124 having a like duration x. Then, cunmt is 
applied in a forward direction with respect to the infiared LED 108 during a third time 
interval 126, also having a duration x. Then, no current is applied to either LED during a 
fourth time interval 128 having a like duration x. Thereafter, the current is again applied in 
the forward direction for the red LED 106 during a fifth time interval 130 which 
corresponds to the first time interval 122. It can be seen that the overall cycle repeats witii 
a period of duration Tequal to 4x. The red LED 106 emits light only when the current is 
z^plied in the forward direction with respect to the red LED 106. Thus, as illustrated by a 
red int^ity waveform 132^ the red LED 106 emits light as a pulse 134 during the first 
time interval 122 and as a pulse 136 during the fifOi time interval 130,andsoorL Thered 
pulses rq>eat with a periodicity equal to T. Similarly, the infi:ared LED 108 emits infiared 
light only when the current is applied in the forward direction with respect to the infiarcd 
LED 108. Thus, as illustrated by an infirared intensity wavefomi 140, the infiared LED 
108 emits infrared light as a pulse 142 during the third interval 126. A next infi:ared pulse 
144 occtirs at an interval T after the infiared pulse 142. Thus, the infiared pulses also 
repeat with a periodicity equal to T. It can be seen that the red pulses and the infrared 
pulses each have a duty cycle of 25%, and the red pulses and the infiared pulses are 
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separated by intervals of one-fourth of each period T the b^innmg of one pulse 
occurs an interval t after the end of the previous pulse). 

As further illustrated in Figure 1, the electromagnetic energy pulses from the red 
LED 106 and the infrared LED 108 are applied to the finger 102. A detector 150 is 
positioned to receive the electromagnetic energy after the energy has passed through a 
portion of the finger 102. The detector 150 is selected to be responsive to both the red 
light and the infrared light and to generate an output signal responsive to the intensity of 
the energy received from each source. An exemplary current ou^ut signal firom the 
detector 150 is represented by a waveform 152 in Figure 2, As illustrated, the detector 
signal waveform 152 conq)rises a first pulse 154 responsive to the first red pulse 134, a 
second pulse 156 responsive to the infrared pulse 142 and a third pulse 158 responsive to 
the second red pulse 136. During the time between the first pulse 154 and the second 
pulse 156, the detector signal waveform 152 comprises noise 160, and during the time 
between the second pulse 156 and the third pulse 158, the detector signal waveform 150 
comprises noise 162. The signal pulses 154, 156 and 158 also include noise superimposed 
thereon. Although shown as repeating noise, it should be understood that the noise varies 
with time. For example, noise caused by ambient light will vary with a periodicity 
corresponding to the 50 Hz or 60 Hz power frequency and their harmonics, particularly 
when the ambient light is provided by fluorescent Ughts which generate significant noise at 
the first harmonic (i.e., 100 Hz or 120 Hz) and the third harmonic (i.e., 200 Hz or 240 Hz). 

The ou^ut of the detector 150 is applied as an input to a signal processor blodc 
170 which processes the detector signal and generates a first signal S^ (/) responsive to the 
detected intensity of the red Ught incident on the detector 150 and generates a second 
signal 5j(0 responsive to the detected intensity of the infrared light incident on the 
detector 150. As illustrated, the signal processing block 170 is synchronized with the LED 
modulator 104 via a set of control lines 180. As wiU be discussed below, the control lines 
180 advantageously communicate signals which provide timing information that 
determines when to activate the red LED 106 and when to activate the infrared LED 108. 
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Figure 3 is a pictorial rqiiesentation of a model of an exemplary system whidi 
incorporates the present invention. The red LED 106 provides a light intensity represented 
as Ird, and the infrared LED 108 provides a light intensity represented as I|r. The e£fects 
of turning the LEDs 106, 108 on and off on periodic bases are modeled by a first multiplier 
or modulator 190 which applies a first modulation signal M\{t) to the red light intensity to 
generate a modulated red signal Irdmod(0 ^ second multiplier or modulator 192 

which applies a second modulation signal Miit) to the infixed light intensity to generate a 
modulated infrared signal Iirmod(0* The modulated light red signal and the modulated 
infrared signal are applied to the finger 102^ or other body portion, as described above. 
The finger 102 has blood flowing therein and is represented in Figure 3 as a block 102. 
The blood in the finger 102 has a volume and scattmng components which vary 
throughout eadi cardiac cycle. The blood carries oxygen and other materials therein. The 
oxygen content is a fimction of both the blood volume and the concentration of the oxygm 
in the blood volume. The concentration of tiie oxygen in the blood volume is generally 
measured as blood oxygen saturation for reasons which are described in fixU in the above- 
identified issued US. Patent Nos. 5,482,036 and 5,490,505. As finther described in the 
two referenced patents, the blood oxygen saturation is determined by comparing the 
relative absoiption of the red light and the infi:ared light in the finger 102. The comparison 
is complicated by the noise caused by movement, ambient light, light scattering, and other 
factors. 

In Figure 3, a pair of signals S\(t) and ^2(0 represent the effect of the.time-varyxng 
volume and scattering conqponents of the blood in the finger 102 on the red light and the 
infi:ared light, req)ectively, passing through flie finger 102 born the LEDs 106, 108 to the 
d^ector 150. The red ligfht signal portion Si(i) is caused by the variable attenuation of &e 
red li^t passing through the finger 102. The infinared light signal portion ^2(0 caused by 
the variable attenuation of the infi:ared light passing through the finger 102. To show the 
effect of the variable attenuations, the signal portion S\iO is illustrated as being applied to a 
first attenuation modulator 191 which multiplies the signal Si(t) by the modulated red 
ou^ut IrdmodW of the first modulator 190. Simflariy, the infi:ared light signal portion Szit) 
is illustrated as being applied to a second attenuation modulator 193 which multiplies the 
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signal 52(0 by the modulated infiared ou^ut Iirmod(0 of the second modulator 192. The 
outputs of the first and second attenuation modulators 191, 193 are provided to the 
receiving photodetector 150. The photodetector 150 is modeled as adder 194 and an adder 
196. The outputs of the first and second attenuation modulators 191, 193 are provided to 
the adder 194 to generate a composite signal A/(/) where: 

M{t) = Sximit) + SiifjMiit). (1) 
The signal A/(0 fro^i the adder 194 is provided to the adder 196 where the signal 
Af(0 is added to a signal n{t) which represents a composite noise signal caused by ambient 
light, electromagnetic pickup, and the like, which are also detected by the photodetector 
150. The output of the adder 196 is a signal Af (0 = + n{f) which includes noise 
components as well as the signal components. The noise components include DC 
components and harmonics of the power line frequency that appear in the ambient light 
In addition, as will be discussed in more detail below, the signal M{t) may also include 
noise at higher frequencies caused, for example, by other devices such as 
electrocauterization equipment, or the like. 

The ^(0 signal ouQ)ut of the third adder 196 (i.e., the output of the detector 150) 
is applied to the input of the signal processing block 170. Withm the signal processing 
block 170, the signal Af (/) is first passed through a fixed gain ampUfier 197 and then 
through an analog bandpass filter 198. The analog bandpass filter 198 has a passband 
selected to pass signals in the range of 20 Hz to 10,000 Hz. Thus, the analog bandpass 
filter 198 removes a significant portion of the noise below 10 Hz. The signal components 
responsive to the blood oxygen saturation are firequcncy shifted by the opwation of the two 
modulation signals M\{() and Mzif) and are passed by the analog bandpass filter 198. 

In the preferred embodiment, tiie output of the analog bandpass filter 198 is 
sampled by an analog-ta<ligital converter 199 and converted therein to digital signals. For 
example, the signals are preferably sampled at 46,875 samples per second. The ou^ut of 
the analog-to-digital converter 199 is a signal MF(k). 

The signal MFik) is provided as a first input to a first demodulating multiplier 210. 
The signal MFQc) is also provided as a first input to a second demodulating multiplier 212. 
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A first demodulating signal D\{k) is provided as a second input to the first demodulating 
multiplier 210, and a second demodulating signal D^ik) is provided as a second input to tfie 
second demodulating multiplier 212. The output of the first demodulating multiplier 210 
is provided as an input to a first lowpass filter 220, and the output of the second 
5 demodulating multiplier is provided as an input to a second lowpass filter 222- The 
bandwidths of the lowpass filters 220, 222 are preferably approximately 10 Hz. 

The output of the first lowpass filter 220 is a signal 5,(0 , which, as discussed 
below, is an estimate of the signal S\{t). The output of the second lowpass filter 222 is a 
signal S^{t) , which, as discussed below, is an estimate of flie signal iSiW- As will be 
10 shown below, the selection of the first demodulating signal D\(k) and the second 
demodulating signal Diik) in accordance with the present invention substantially reduces 

or eliminates the effects of noise in the two output signals and S^{t) and also 
substantially reduces or eliminates crosstalk between the two signals. 

In the preferred embodiment of the present invention, the sample rates of the 

15 outputs of the lowpass filter 220 and the lowpass filter 222 are compressed by respective 
sample rate compressors 221 and 223. In particular, the sample rate compressors 221, 223 
reduce the sample rate by 750 to a sample rate of, for example, 62.5 Hz to provide an 
output which can be further processed in accordance with the methods and apparatuses 
described in the above-referenced patents. The sample rate compressions which occur in 

20 the sample rate compressors 221, 223 reduce the rate at which the output signals iS, (0 and 
^2^) need to be processed while maintaining the sample rate well above the 0-10 Hz 
frequency content of the signals of interest The outputs of the filters 220, 222, or the 
sample rate compressors 221, 223, if included, are provided on respective output lines 224 
and 226. 

25 In order to facilitate an understanding of how the present invention operates in 

demodulating the ou^ut signal MF{k) torn the analog-to-digital converter 199, the 
modulation signals M\{t) and Miit) will first be described in terms of their frequency 
components. One skilled in the art will appreciate that the modulation signals M\it) and 
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A/2(0 <^ ^® represented as a Fourier cosine series expansion 

(^•8->5]*^ a^cos{n<ot), where g> = In/T) representing the fundamental and haimonic 
fiequencies of the rectangular signal pulses. One skilled in the art will luiderstand that the 
Fourier scries expansion includes phases; however, by suitably selecting the time origin, 
the phases are set to zero- A component which is ISO*' out of phase with a corresponding 
component will advantageously be represented by a minus sign before the coefficient 

Figure 4 illustrates a frequency spectra of the first modulation signal M\{t) for n « 
0, 1, 2, . . where the horizontal axis represents frequency, with the energy in the DC 
component along the vertical axis and increasing hamionics of the fundamental fi^ucncy 
along the horizontal axis. The length of each component ofM\{() along the vertical axis 
represents the energy E(n) in each component of the frequency spectra. The first 
component to the right of the vertical axis is at the fundamental fi^uency (i.e., l/T), which 
is designated herein as fo; however, it should be understood that the fundamental frequency 
fo corresponds to n = 1. The second component to the right of the vertical axis is the first 
harmonic/i (i.e., n = 2), which has a fi:equency which is twice the fimdamental firequency. 
The third component to the right of the vertical axis is the second haimonicTi (i.e^ n =3), 
which has a frequency which is three times the fimdamental frequency. The components 
to the right of the second harmonic are nimibered accordingly. (Note, other conventions 
identify the fundamental frequency as the first harmonic, and designate the second 
harmonic as the frequency that is twice the fimdamental frequency. The identification of 
the fundamental firequency as^o is used in the discussion that follows.) 

In Figure 4, a modulation envelope 230 is shown in dashed lines. The modulation 
envelope 230 represents the magnitudes of the fimdamental and the harmonics of the 
signal M\(i). the shape of the envelope is determined by the modulation signal M\(t) 
which, for a repeating rectangular pulse train starting at time / - 0 and having a normalized 
amplitude of 1 , can be expressed as: 
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Where sine is the function (sin 70c)/7ix (i.e., sinc(nT/7) = sin(mrc/I)/(mrc/7)). In the 
example shown, x = */ir. (Note that for sampled signals, the envelope is more accurately 
represented as sino/sinP; however, as well known in the art, for the frequencies of int^est, 
the sine function is a suitable approximation.) Thus, the frequency spectra has nulls at 
5 n = 4, n = 8, /i = 12, and so on, corresponding to the third harmonic yi, the sevoith 
harmonic /?» the eleventh harmonic /u, and so on. Note that Equation 2 is an idealized 
form of the equation for M\(t% and that in general: 

Mi(t)-^±a,sinc^e'^^' (3) 

where a„ is a complex numbo*. In the discussion that follows, the values of an are assumed 
10 to be real numbers only. 

A similar frequmcy spectra (not shown) for the modulation signal M2(t) is 
d^ermined by the expression: 

M,ry»fgr-//sinc[^]cos[i^] (4) 

An envelope for the frequency spectra of second modulation signal M2it) will have 
15 the same magnitudes; however, it should be understood that because of the (-1)" temi in 
the expression for Miit), the fundamental fo and every even harmonic (i.e.,yi,y4, etc.) are 
180"^ out of phase with the corresponding harmonic of the first modulation signal M\(t). 

In Figure 3, the analog-to-di^tal converter 199 converts the signal Af(t) to a 
sequence of sampled digital values MF{k) at a sampling rate of; for sample, 46,875 
20 samples per second. As discussed above, the first demodulating niultiplier 210 multiplies 
the output MFlk) of the convater 199 by the first demodulating signal D\{k) to generate 

the first output Sequence {k) , and the second demodulating multiplier 212 multiplies the 
output MF(k) by the second demodulating signal Diik) to generate the second output 

sequence Sj (k) . The multiplication by the multipliers 210, 212 can also be expressed as 
25 follows: 

S,ik)=LP[MF{k)Di{k)] (5) 

and 

.16- 



Substitute Sheet (Rule 26) 



o o 

wo 99/52420 PCTAJS99/07825 



S,(k) =LP[MFik)D2ik)] (6) 
where LP is the transfer function of the lowpass filter 220 and of the lowpass filter 222. 
If, for simplicity, the noise is assumed to be zero, then: 

Af (0 = SiimiO + S2(/)M2(0 (7) 

ThCTcfore: 

5, (k) = LP[ [SiikWiik) + S2ik)M2ik)]Diik) ] (8) 

and thus 

5, (A) - LP[ [SiikWximDiik) + lS2ik)M2{k)]Dm ] (9) 

Similarly: 

5,(Jfc) = LP[ [S2(k)M2ik)]Di(k) + [SxikWmVhik) ] (10) 
Since LP is a linear operator, the right-hand side of Equations 9 and 10 can be 
split into two terms. The first terai on the right-hand side of each of Equations 9 and 10 
above is the desired signal portion of the equation, and the second terai on the right- 
hand side of each of the equations is the crosstalk portion. Thus, in order to reduce the 
crosstalk to zero, the second term of each of Equations 9 and 10 is set to zero: 

LP[ S2{k)M2ik)Dx(k) ] = 0 (11) 

and 

LP[ SxikWiikWk) ] = 0 (12) 
By setting the second terais to zero. Equations 9 and 10 reduce to: 

S, (it) = LP[ Sx(k)Mxik)Dx(k) ] (13) 

and 

S,ik)^LP[S2{k)Miik)D2(k)] (14) 
One goal of the present invention is to select the demodulating signals Dx{k) and 
D2(k) to satisfy Equations 1 1 and 12 to thereby reduce Equations 9 and 10 to Equations 13 
and 14. This is accomplished by utilizing Equations 2 and 3 to simplify the two equations 
by selectively using components of the two modulating signals Mx{t) and A/2(0 ^ generate 
the demodulating sequences Dxik) and Djik). 

In order to simplify the discussion. Equation 2 can be rewritten as: 
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Mi(t) - 2;£r«;cosr/ia)r; (is) 

where E(n) is the sine envelope for the fundamental frequency To = 1) and the haimonics 
fx {n = 2Xf2 (n - 3), and so on, where cos(n(ot) represents the cosine term cos(27rn//7), 
where co=27i/r. (Note, as discussed above, for discrete sampled signals, the actual 
envelope of £(n) is a sina/sinp function; however, for the frequencies of interest, the sine 
5 function is a suitable representatioa) 

As discussed above, the DC term (n - 0) does not need to be considered because of 
the operation of the filter 198, and the anaiog-to-digital converter 199, as well as the action 
of the demodulation, which shift any unwanted DC or low frequency signals having a 
frequency less than approximately 10 Hz (hereinafter near-DC signals) to higher 
10 frequencies before lowpass filtering. As a further simplification, the magnitude of tfie 
fundamental term in Equation 15 is normalized to a value of 1 (i.e., E(l) = 1). Note that 
the normalization results in the need for a scale factor, which will be discussed below. 
Thus, Equation IS becomes: 

A/|(0 = cos ©r + a cos 2(ot + b cos 3a>t + c cos 4(ot + • • • (16) 
15 The demodulation signal D|(0 is defined as: 

Z)|(0 = cos ©/ + B cos 2fi)r (17) 
For reasons set forth below, only the first two cosine terms are needed. 
Similarly, the second modulating signal Mtit) becomes: 

A/2(0 - -cos ©r + a cos lot - b cos 3©/ + c cos 4©/ + . . . (18) 
20 and the second demodulating signal I>2(0 is defined as: 

Z)2(0 = -cos©/ + Bcos2©r (19) 
Note that the signs of the fundamental and odd hamionics in Equation 18 are ISC' out of 
phase with the corresponding terms in Equation 16. 

Note, as will be developed more fully below, by including only the fimdamratal 
25 (cos ©r) and the first harmonic (cos 2©r) in each of the demodulation signals, only tfie 
signals proximate to the fundamental and first harmonic need to be considered. By 
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eliminating higher harmonics, the effects of the higher harmonics of the power line 
frequency are also eliminated in the output signals generated by the present invention. 

Assume that the filter 198 and the analog-to-digital converter 199 do not affect the 
magnitude of the signal MFik) with respect to Af(t) for the fiequendes having significant 
energy. Therefore, starting with Equation 7 above, Af(t) can be written as: 
Af(t) = iSi(0[cos cor + a cos 2©/ + b cos 3(or + . • •] 

+ SiiOl-cos ©/ + a cosZco/ - b cos 3©/ + . . .] (20) 
When the first demodulating multiplier 210 multiplies Afit) by D\(t), the terais on the 
right-hand side of Equation 20 are multiplied by the terms on the right-hand side of 
Equation 17. Thus: 

Af (r)Di(0 S\ (/)[cos (0/ + a cos 2©/ + b cos 3©/ + . . .][cos ©f + B cos 2©/] 

+ SiiOl-cos ©/ + a cos 2©r -b cos 3©r + . . .][cos ©/ + B cos 2©r] (21) 
The term 5i(r)[cos ©/ + a cos 2©r + b cos 3©r + . . .][cos ©/ + B cos 2©/] is the signal temi 
which is to be preserved, and the terai SiiOi-cos ©r + a cos 2©r -b cos 3©/ + . . .][cos ©r+ 
B cos 2©r] is the crosstalk temi to be eliminated 

Expanding the crosstalk terai from Equation 21, generates: 

crosstalk = SiiOi-cos^m-B cos ©f cos 2©/ + a cos2©/ cos ©r + aB cos^©r 

- b cos 3©r cos ©r - bB cos 3<ot cos 2©r + . . .] (22) 
Using the identity, cos(x)cos(y) = l4[cos(x+y)+cos(x-y)], the crosstalk term 

from Equation 22 becomes: 

crosstalk = S2{i)[ -Yiicos la^t + I) + ((a - B)/2)[cos 3©r cos ©r] 
+ (aB/2)[cos 4©/ + 1] - (b/2)[cos 4©f + cos 2©f] 

- (bBy2)[cos 5©r + COS ©/] + .• .] (23) 
The remaining Xeax& in Equation 23 will all have a factor of cos ©/ or higher. Thus, 
Equation 23, when fiilly expanded only includes near-DC terms: 

crosstalkDc = LP[S2(0[(aB/2)-*/2]] (24) 
where ^2(0 conresponds to the infirared portion of the original plethysmogr^h signal 
which has a bandwidth of interest of approximately 0 to 10 Hz. Any components present 
above 10 Hz will be eliminated by the action of the lowpass filter 220. Thus, it can be 
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seen that only the signals of interest are folded back to DC or near-DC. By using the 
lowpass filter 220, the DC teims and near-DC terms can be isolated so that only the DC 
terms and near-DC terms of the crosstalk are presented at the output of the lowpass filter 
220. Thus, in order to eliminate the crosstalk, the crosstalk terms in Equation 24 need to 
be set to zero: 

LP[52(0[aB/2-/2]]=0 (25) 

Thus: 

B = l/a (26) 
The result in Equation 26 can alsb be expressed using a geometric interpretation of 
vector projection (i.e., dot products) of ^2(0 and S\(t) wherein the projection of ^2(0 onto 
D\(t) is equal to zero and the projection of -S2(/) onto DiQ) is maximized. In other words, 
express iSt(0> Szit^ D\(t) and Dzii) as vectors of samples in an n-dimensional sample space 
(e.g.. Slit) is represented as a vector Si of samples Si{k)). For example, in a prefened 
embodiment, n = 148, and thus Sj, S2, Di and D2 are vectors of 148 samples each. The 
first crosstalk term is SrD2. The second crosstalk term is S2*Di. The first signal output is 
Si-Di. The second signal output is 82-02. Select the vectors Di and D2 to drive the 
crosstalk toms to zero. 

The relationship in Equation 26 also works to preserve the signal temi. In 
particular, the signal term in Equation 21 can be expanded and lowpass filtered in the same 
manner as the crosstalk term to obtain: 

signal^ 5,(/) = LP[Si(0[(aB/2) + '/2]] (27) 
Using the relationship fiom Equation 26, then Equation 27 becomes: 

signal = S, (0 = LP[5i(/)[(a/2a) + Vi] = LP[5i(r)] = S,(0 (28) 
It can be readily shown that the same relationship holds for the crosstalk term and 
the signal term for the signal ^2(0 by defining the second demodulation signal Diit) as: 

jD2(/) = -cos<Dr + Bcos2(0^ (29) 
and multiplying MiQ) by Dtit). After expanding the crosstalk and signal tenns and 
eliminating the trans above 10 Hz, it can be diown that by selecting B- l/a, the crosstalk 
term is canceled and the signal term ^2(0 is recovered. 
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From the foregoing, it can be seen that by choosing the relationship between the 
magnitude of B as the reciprocal of a, then the crosstalk terms are eliminated and the signal 
temis are preserved. Note that neither A nor B is an absolute value. As set forth in 
Equation 16, a is the magnitude of the cos2cor term of M\{t) when the magnitude of the 
cosco/ term of A/|(0 is nomialized to 1. Similarly, from Equation 17, B is the magnitude of 
the cos 2(or term of Z)i(0 when the cos cor term o{D\(t) is normalized to 1 . 

It should be understood that both D\{t) and DiQ) can include higher harmonic 
terms; however, such additional teraas could result in increased sensitivity to the noise of 
fluorescent lights and the like because of the haraionics of the 60 Hz power line frequency 
(or the 50 Hz power line frequency in other countries). For ©cample. Figure 5 illustrates 
an exemplary spectrum of the first and second harmonics of the present invention when the 
fundamental fi^uency is selected to be 3 16,7 Hz. Thus, the first harmonic frequency is 
633.4 Hz. Note that the variations in the signals caused by blood flow throughout a 
cardiac cycle causes the fundamental and harmonics of modulation frequency to be 
surrounded by sidebands representing the frequency content of the plethysmograph. For 
example, in Figure 5, the first and second harmonics are at 316.7 Hz and 633.4 Hz, 
±10 Hz. 

As further illustrated in Figure 5, the conventional 60 Hz power line fi^uency has 
harmonics at 120 Hz, 180 Hz, 240, etc. Thus, the nearest harmonics of the power line 
frequency to the first harmonic of the present invention are at 300 Hz and 360 Hz, and the 
nearest harmonics of the power line frequency to the second harmonic of the presmt 
invention are 600 Hz and 660 Hz. Similaily, if used in a country having a 50 Hz power 
line frequency, the nearest harmonics to the first harmonic of the present invention arc 300 
Hz and 350 IJz, and the nearest harmonics to the second harmonic of the present invention 
are 600 Hz and 650 Hz. Evenif thepowerfiiequcncy wereto varyby up to 1.5percOTt, 
the noise generated by the ambient light finom fluorescent lamps, or the like, would not be 
at the firet and second harmonic frequencies of the present invention. The fundamaital 
fi^uency has thus been selected to avoid power line caused ambient noise at the first and 
second harmonic fi:equencies. 
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The foregoing discussion assumed that the filter 1 98 did not sigbificantly affect the 
amplitude of the filtered signal. If the filter 198 does have an affect on the amplitude, then 
B will be a constant times the value of B detemiined above: 

B = k/a (30) 
5 where k depends on the relative attenuation of the first hannonic and the second harmonic 
throughthe filter 198. 

Although the value of the coefficient B can be calculated as set forth above, the 
calculations may be complicated if the filter 198 or the modulators 190, 192 introduce 
phase changes which cause the calculatibns to be performed on complex numbers* For 
10 exanq)le, if the modulation signals M\{t) and Mtit) arc not rectangular waves wWch have 
25% duty cycles and which are precisely 180'' out of phase, as illustrated herein, thm the 
coefficients of the fiequency components of the modulation signals may be complex to 
accoxmt for the phase relationships, and thus, the coefficients of the demodulation signals 
may be complex. 

15 As illustrated in Figure 6, the value of B can also be determined empirically by 

performing a initial measurement with one channel (i.e., either the red pulse or the infirared 
pulse turned off) and minimizing the crosstalk. In particular, during the initial 
measurement, the waveform 140 in Figure 2 is set to a continuous zero value so that no 
infiared pulses are generated. Thus, tiie detector 150 (Figure 1) receives only tiie light 

20 generated by the red LED 106. Thus, MiQ) is set to zero, and Equation 10 for 5^(0 
becomes: 

5,(r) = LP[ Sximmiit) ] (31) 

It can be seen that SjCO includes only a crosstalk portion, which can be measured on the 
output from tiie second lowpass filter 222. Thus, by varying tiie value B while monitoring 

25 the magnitude or tiie RMS (root-mean-squared) value of the output signal S^iO^ a 
minimum magnitude (0 mm for tiie ou^ut signal (0 can be found which corresponds 
to tiie best value Bbest for B . In an ideal system, the best value for B corresponds to a zero 
value for the output signal (r) ; however, in a real environment, tiie best value of B may 
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correspond to a non-zoo value for 5, (0 (i-e., a minimum error for 5, (0 ). It should be 
understood that the value of Bbest can also be determined by turning off the red LED 106 
and varying B while monitoring 5,(0 until 5,(0 is minimized. 

From the foregoing, it can be seen that the effect of the modulation signals D\(t) 
and DiiO is to shift the DC or near-DC noise terms up in frequency while shifting the 
signals of interest at the harmonics back to DC or near-DC, which in effect interchanges 
the noise spectra and the signal spectra so that the noise spectra can be eliminated by the 
action of the lowpass filters 220, 222, leaving only the signals of interest 

Figure 7 illustrates a prefored embodiment of the present invention wfaidi 
implements the fimctions described above in a digital system. Prefoably, the digital 
system comprises a digital signal processor (not shown), and the blocks d«cribed herdn 
comprise data structures within the digital signal processor and software routines that 
implement the processes described below. In particular, the present invention comprises 
an LED demodulation block 300 which receives a digital configuration signal on a bus 
310, a clock signal on a line 312 and a digital detector signal on a bus 314 as inputs. The 
digital configuration signal bus 310 provides a way to change the configuration of the LED 
demodulation block 300 to accommodate different LEDs and different detection 
algorithms. Preferably, the clock signal on the line 312 is a 46,875 Hz (46.875 kHz) , 
square wave signal which is used to synchronize the timing fimctions of the present 
invention. The digital detector signal on the line 314 is the output of the analog-to-digital 
converter 199. The analog-to-digital converter 199 is cormected to the ou^ut of the 
detector 150 (via the amplifier 197 and tiie filter 198) and samples the output of the 
detector 150 at 46,875 samples per second to provide a stream of sampled digital values of 
the red light arid infiared light incident on the detector 1 50. 

The LED modulation block 300 generates a demodulated red signal output on a 
bus 340 and generates a demodulated infirared signal output on a bus 342. TTie 
demodulated red signal output is passed through the low pass filter 220 and is output 
therefrom as the signal S,(0 • "Hie demodulated infiared signal output is passed through 
the low pass filter 222 and is output tiierefiom as the signal ^,(0 • As fiuther illustrated in 
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Figure 8, the LED demodulation block 300 comprises a modulo-M block 3S0» an LED 
demodulation state table block 352, the first demodulating multiplier 210 and the second 
dCTiodulating multiplier 212* 

The modulo-M block 350 receives the main 46,875 Hz clock signal on the tine 312 
5 as one input and receives a MODULUS signal on a bus 354 as a second input The bns 
354 fomis a portion of the configuration bus 310. The modulo-M block 350 divides the 
clock signal by the MODULUS signal and generates a RESIDUE signal (described below) 
on a bus 356 which is provided as one input to the LED modulation state table block 352. 
The LED modulation state table block 352 also receives the configuration signals on the 
10 configuration bus 310. 

The LED demodulation state table is responsive to the residue signal and the 
configuration signals to generate the first demodulating signal D\it) on a bus 360 and to 
generate the second draiodulating signal Diit) on a bus 362. The first demodulating signal 
jD|(0 is provided as one input to the first demodulating multiplier 210, as described above. 
15 The second demodulating signal Jhif) is provided as one input to the second demodulating 
multiplier 212, as described above. The first demodulating multiplier 210 and the second 
demodulating multipUer 212 receive the digital detector signal on the line 314 as 
respective second inputs. The demodulating multipUers 210, 212 multiply the digital 
detector signal by the first demodulating signal D\{t) and the second demodulating signal 
20 DiiOj respectively, to generate a demodulated red signal and a demodulated infixed signal 
on the buses 340 and 342, respectively. Because the outputs of the two demodulating 
multipliers 210 and 212 include the temos cos co/, cos 2g)/, and higher, the demodulated 
signals on the buses 340 and 342 are provided as respective inputs to the low pass filters 
220 and 222 to pass only the near-DC terms, as discussed above. The outputs of the 
25 lowpass filters 220 and 222 on the buses 344 and 346, respectively, are the 5, (/) signal 

and the (/) signal which contain only the near-DC terms, vAActu in accordance witii the 
discussion presented above represent the ori^nal input signals S\(t) and ^2(0 with the 
unwanted noise substantiaUy reduced or eliminated. The two signals 5,(0 and S^it) are 
then ^lied to computation ciicuitiy (not shown) which computes the blood oxy^ 
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saturation and other cardiographic parameters in a manner described in the above-cited 
U,S. Patent Nos, 5,482,036 and 5,490^05, 

The residue signal generated as the output from the modulo-M block 350 is a 
multiple bit signal that counts from 0 to MODULUS-1. In the preferred embodiment 
described herein, MODULUS has a value of 148. Thus, the RESIDUE output of the 
modulo-M block 350 counts from 0 to 147. The RESIDUE output of the modulo-M block 
350 is a number that is provided as the input to the LED demodulation state table block 
352, As illustrated in Figure 9, the RESIDUE output on the bus 356 corresponds to the 
signal 180 in Figure 1 and is also provided to the input of an LED modulation state table 
block 370 which, together with an LED driver circuit 372, comprise the modulation block 
104 (Figure 1) which generates the drive signals to the red LED 106 and the infrared LED 
108. As described above, the red LED 106 and the infrared LED 108 genwate the 
modulation signals M\{t) and Mi{i), respectively, which effectively operate as carriers for 
the plethysmograph waveform to be measured. In particular, as illustrated by a red drive 
timing waveform 374 and by a infrared drive timing waveform 376 in Figure 10, the 
modulation state table block 370 generates a red signal pulse 378 during the time when the 
RESIDUE signal increments from 0 to 36. Then, tfie modulation state table block 370 
generates neither a red signal pulse nor an infrared signal pulse during the time when the 
RESIDUE signal increments from 37 to 73. Then, the modulation state table block 370 
generates the infrared signal pulse 380 during tiie time when the RESIDUE signal 
increments from 74 to 1 10. Then, the modulation state table block 370 again generates 
neither a red signal pulse nor an infrared signal pulse during tiie time when the RESIDUE 
signalincrCTicntsfiomlll tol47. The RESIDUE signal then resets to 0 and the process 

repeats continuously. 

The red signal pulse 378 and the infrared signal pulse 380 from the modulation 
state table block 370 are provided as inputs to the LED driver circuit 372 which turns on 
the red LED 106 when the red signal pulse 376 is active and turns on the infrared LED 108 
when the infrared signal pulse 378 is active by gmerating the current waveform 120 
illustrated in Figure 2. The circuitry for convating tiie red signal pulse 376 and the 
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infiared signal pulse 378 to the bi-directional current pulses of the waveform 120 is 
conventional and does not need to be described herein. 

In the preferred embodiment, the LED demodulation state table block 352 
implements demodulation equations which generally correspond to the Equations 17 and 
19 described above. In particular, the LED demodulation state table block 352 receives 
the RESIDUE as one input to the state table and steps through the state table based upon 
the current value of the RESIDUE. The LED demodulation state table block 352 
generates two output values for each value of the RESIDUE, wherein the first output value 
is the first demodulation signal A(0 the signal bus 360, and the second ou^ut value is 
the second demodulation signal Di{f) on the signal bus 362. 

In particular, the LED demodulation state table block 352 implements the 
following forms of the demodulation signal D\{t) and the i>20D equations: 

/ r (32) 

and 

In Equations 32 and 33, the value SCL is a scale factor which detennines the magnitudes 
of the two demoduladon signals and which is used to compensate for the normalization 
discussed above and to compensate for other factors, such as, for example, non-ideal 
rectangular pulses. The method ofdeterniining the scale factor will be set forth below. In 
one particularly preferred embodiment, the value of SCL is 2221441469. The value 
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HWD is a hardware distortion factor, whidi corresponds to the value of B discussed 
above. The determination of the value B was described above, and will be described again 
below in connection with this prefenred embodiment. In one particulariy preferred 
embodiment where the pulses applied to the red LED 106 and the infiared LED 108 are 
idealized rectangular waves having 25% duty cycles, the value of HWD can be calculated 
to be 1.414213562. This ideal value for HWD can be determined by recognizing that the 
value of the coefficient A for the cos 2(ot tenns in Equations 16 and 18 is determined by 
the sine fiinction. When the coefficient of the cos of temi is nomuilized to 1 , as in tiie two 
equations^ then the value of the coefficient a is equal to VJ/2 • Hius, the ideal value for B 
(i.e., HWD) is V2. Of course, the actual value of the coefficient B, and thus HWD, will 
vary when the red pulses and the infrared pulses are not true rectangular waves. Since, in 
actual embodiments, the pulses will have finite rise times and fall times, the optimum 
value of HWD is preferably found empirically in the manner described below. 

The value 18.5 in Equations 32 and 33 is used to align the demodulation 
waveforms with the modulation waveforms so that the peak of the cosine fimctions 
corresponds to the midpoints of each of the modulation waveforms. The value HWA is a 
hardware delay factor which may be needed in certain embodiments to compensate for 
delays in the analog processing, the digital processing or both, which cause the 
demodulation signals D\(t) and 2)2(0 to be out of phase with the modulation signals Mt(t) 
and Miit). In an ideal environment, flie value of the hardware delay fector is 0. However, 
in one particularly preferred embodiment, tfie value of the hardware delay factor is 39. 
The modulus was described above and is basically the number of steps in each period of 
the waveforms. In the embodiment described herein, the modulus is 148. The value R is 
the residue! which varies from 0 to modulus-1, and thus, in the preferred embodiment, 
R varies from 0 to 147. 

In operation, the clock signal on the line 312 causes the modulo-M block 350 to 
generate the RESIDUE signal, as described above. The RESIDUE value is applied to the 
LED modulation blodc 104 which generates the modulation signals Mi(0 and Mjit), as 
described above. The RESIDUE value is also applied to the LED demodulation state table 
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block 352 which generates a new value for D\it) and a new value for Diit) for each new 
RESIDUE value. Thus, 148 values o{D\(t) and DzQ) are generated for each complete 
cycle. Because the clock signal is operating at 46,875 Hz, the modulation signals M\{t) 
and MiiO and the demodulation signals D\(0 and Diit) have a fimdamental frequency of 
5 316.722973 Hz, which, as discussed above, does not correspond to any harmonic of 
conventional 50 Hz or 60 Hz powa: line frequencies. 

The HWA (hardware delay factor) value, the HWD (hardware distortion factor) 
value and the SCL (scaling factor) value are found empirically as follows. First, the ideal 
values of the hardware delay &ctor, the hardware distortion factor and the scale factor are 
10 applied to the Equations 32 and 33 in the LED demodulation state table block 352 (i.e., 
HWA = 0, HWD = 1.414213562, and SCL = 2.221441469). To determine the optimum 
value of the hardware delay factor, the second modulation signal MiQ) is set to a constant 
value of zero (i.e., the infrared LED is maintained in its OFF state). The red LED pulses 
are applied as set forth above, and the digital detector output signal from the analog-to- 
15 digital converter is monitored and compared to the modulation signal M\{t). The relative 
delay between the beginning of the modulation signal M\(t) and the detection of the 
beginning of the responsive output from the analog-to-digital converter is the optimum 
hardware delay factor (HWA) value. In one exemplary embodiment, the optimum value of 
the hardware delay factor is 39. 
20 After determining the value of the hardware delay factor and applying it to 

Equations 32 and 33, the ideal value of the hardware distortion factor and the ideal vahie 
of the scale factor are applied to the two equations. Again, with the red LED pulses 
applied to the red LED 106 and no pulses applied to the infrared LED, the value of the 
hardware distc^on factor is slowly varied from its ideal value v/bHc the DC component of 
25 tiie demodulated infrared signal output on the line 342 is monitored. The value of the 
hardware distortion factor is varied until the measured DC component is minimized, and 
the value of the hardware distortion factor corresponding to the miiumal DC component is 
selected as the optimum value for the hardware distortion factor. 
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Next, with the value of the hardware delay factor and the value of the hardware 
distortion factor set to their respective optimum values, as determined above, the value of 
the scale factor (SCL) is initially set to 1. Again, with the modulation system generating 
pulses only to die red LED 106, the DC component of tfie demodulated red signal output 
on the line 340 is measured. In addition, the dififercnce in amplitude between the on state 
and the off state of the digital detector signal from the filter 198 is measured. The ratio of 
the measured amplitude difference to the measured DC component of the demodulated red 
signal ou^ut is selected as the optimum value for the scale factor. 

An exemplary demodulation waveform Dx(t) is illustrated by a waveform 400 in 
Figure 11 and an exemplary danodulation waveform Diit) is illustrated by a waveform 
402 in Figure 11. The demodulation waveforms in Figure 11 are illustrated with the 
hardware delay factor set to 0 in order to align the waveforms with the modulation 
waveforms in Figure 10. It should be understood that when the hardware delay factor is 
non-zero, the demodulation waveforms in Figure 1 1 will be shifted in phase with respect to 
the modulation waveforms in Figure 10. 

Although described above in coimection with the variation of the amplitude of the 
first harmonic component of the demodulation signals in order to minimize the crosstalk, it 
should be understood that the relative amplitude of the second harmonic component of the 
demodulation signals witii respect to the amplitude of the fundamental component of the 
demodulation signals is detemiined by the relationship of the amplitude of the first 
harmonic component of the modulation signals to the ampUtude of the fimdamental 
component of the modulation signals. Tne relationship of the amplitude of the first 
harmonic conqwnent of the modulation signals depends in part upon the duty cycles of the 
modulation signals^ If the modulation duty cycles are varied, the amplitude of the first 
harmonic component of the modulation signals changes. Thus, the crosstalk may also be 
minimized by holding the amplitudes of the components of the demodulation signals 
constant while varying tfie duty cycles of the modulation signals. One skilled in the art 
will appreciate that other variations in the modulation and demodulation signals may also 
be used to minimize the crosstalk between the two output signals. 
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A plurality of signals Su Si^ . . . Sn canhe demodulated and the crosstalk 
between signals reduced to a minimum by application of the foregoing invention to more 
than two signals. 

Additional information can advantageously be derived from the digitized detection 
signal on the bus 314 and can be used to provide indications regarding the reliability of the 
demodulated signals generated as described above. In particular, although the present 

system is capable of demodulating the iS, (t) signal and the S2 (t) signal in the presence of 
significant ambient noise from light and other sources, it is possible that the level of tiie 
ambient noise is sufficiently high to affect the demodulated signals. Figures 12 and 13 
illustrate a time domain method and system for detennining the ambient noise level, and 
Figures 14 and IS illustrate a frequency domain method and system for determining the 
ambient noise level. 

As illustrated in Figure 12, the digital detection signal 152 is sampled by a sample 
signal represented by a waveform 500, which comprises a plurality of sampling pulses 
502. The sampling pulses 502 are timed to occur during the intervals between the red 
pulses 134, 136 and the infrared pulses 142, 144 when no red light and no infrared li^t 
should be detected by the detector 150 (Figure 1). Thus, any energy detected during the 
sample intervals is primarily caused by ambient light and other noise soiuces. As 
illustrated, the sampling pulses 502 preferably occur at the approximate mi(^oint of each 
interval between the red and infrared piilses. 

As illustrated in Figure 13, the digital d^ection signal bus 314 is provided as an 
input to a time doinain sample 520. The time domain sampler 520 also receives tiie 
RESIDUE signal on the bus 356 as a second input Tlie time domain sanq>ler is 
responsive to the RESIDUE signal to sample the digital detection signal at times when tiie 
value of die RESIDUE signal corresponds to the quiescent times of the red pulses 134, 136 
and the infrared pulses 142, 144, As described above, the red pulses 134, 136 are 
generated when the RESIDUE signal has values between 0 and 36, and the infrared pulses 
are generated when the RESIDUE signal has values between 74 and 1 10. Thus, assuming 
no hardware delay, the sampling pulses 502 are preferably generated, for example, when 
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the RESIDUE signal has a value of 55 and when the RESIDUE signal has a \-aIue of 129, 
which positions the sampling pulses at the approximate midpoints of the quiescent 
intervals between the pulses. As discussed above, the actual system has a hardware delay 
caused by processing times. Thus, if the system has a hardware delay factor of, for 
example, 39, the sampling pulses 502 are shifted in time to occur when the RESIDUE 
signal has a value of 94 and a value of 20 (168 moduloi48). The sample times used by the 
time domain sampler 520 are advantageously determined by configuration signals received 
via the digital configuration bus 310, described above. For example, the time domain 
sampler 520 can be initiaUy set to sample at RESIDUE signal values of 55 and 129, and 
the value of the hardware delay value factor (HWA) communicated by the digital 
configuration bus 310 is added to both values to shift the sample to the correct sample 
interval. 

As illustrated in Figure 14, a detection signal spectra 550 includes the two 
fi-equency components corresponding to the flmdamental and first harmonic of the 
modulation signal at 316.7 Hz and 633.4 Hz, respectively. The spectra 550 further 
includes the fimdamental and multiple harmonics of the 60 Hz power line fiequency. In 
addition, tiie spectra 550 includes noise at a raultitiide of frequencies which may be caused 
by various sources. One particularly tit>ublesome source of noise encountered in pulse 
oximrtry systems is an electrocauterization device, which uses a high frequency electrical 
current to make surgical incisions and to cauterize the surrounding blood vessels at the 
same time. Although primarily high fi^juency noise sources, such devices also generate 
significant noise at lower frequencies because of arcing. When an eledrocauterizarion. 
device is operated close to a pulse oximeter detector, the noise generated by the device can 
overwhelm the signals generated by the pulse oximetry detector. In other words, the noise 
floor can be greater tiian the detectable signal fi»m the pulse oximetry detector. 

It is desirable to detect when the noise floor is too high so that the pulse oximetry 
system can indicate that the demodulated signals may not be reliable. In order to 
determine the level of the noise floor, the present invention samples the spectra 550 to 
detennine the content of the frequency components deterted at frequencies other tiian the 
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fundammtal and harmonic frequencies of fhe modulation signals. In paiticular, as 
illustrated by a sample control signal 560 in Figure 14, the portions of the spectra 550 
which do not include the fundamental and harmonics of the modulation signal are 
sampled. Thus, in the preferred embodiment, the magnitudes of the spectra at 316.7 Hz, 
633.4 Hz, 950.1 Hz, etc., are not sampled. Furfliermore, because a band of frequencies 
around the fundamental and harmonics of the modulation signal also include significant 
infomiation caused by the modulation of the red pulses and the infrared pulses by the 
changes in blood flow during each cardiac cycle. Thus, as illustrated in Figure 14, in the 
preferred embodiment, a band of frequencies surrounding the fundamental and harmonic 
frequencies of the modulation signals (i.e., the sidebaruls discussed above) are not included 
in the samples. For example, a band of at least ±10 Hz around each of the fundamental 
and harmonic frequencies is not included in the sanq>les. 

The intensities at the sampled frequencies are averaged, and an ou^ut signal is 
generated which represents the average inte^ty of the noise signals. Other portions (not 
shown) of the digital processing system advantageously monitor the average intensity of 
the noise signals, and, if the average intensity exceeds a selected threshold based upon the 
size of the measured plethysmograph, then the demodulated output signals from the 
system are considered as being unreliable and shoxild not be used. 

Figure 15 illustrates a preferred embodiment of a system that determines the noise 
floor, as described above. The system of Figure 15 includes a Fast Fourier Transform 
block 600 which receives a plurality of samples from the digitized detector bus 314 and 
generates a transformed output on a bus 610. The transformed output on the bus 610 
rqiresents the qiectra of the samples. In the preferred embodiment, a sufficient number of 
samples are taken to rq^resent approximately 44 milliseconds of data so tfiat at least two . 
cycles of the 60 Hz power are included withm the samples. For example, appioxlmatdy 
1,024 samples can be taken during the 44-millisecond interval at a sanq)Ie rate of 
approximately 23.4 kHz (e.g., one-half the systrai timing rate). The spectra for a 44- 
millisecond interval are provided as inputs to a spectral sampler 620 ^ch eliminates- the 
samples in the ^10 Hz bands around the fundamental and harmonic frequencies of flie 
modulation signals. The output of the spectral sampler 620 is provided on a bus 630 and is 
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thereby provided as an input to an averager 640. The averager 640 averages the sampled 
noise spectra which it receives and provides an averaged output on a bus 650. The 
averaged output on the bus 650 represents the noise floor and is provided to other portions 
of the digital processing system where it is compared to the selected threshold to determine 
whether the noise floor is excessive. The threshold is not necessarily fixed, but is 
dependent on the strength of the plethysmograph, which in turn depends upon the 
perfusion of blood in the body portion being measured. 

The embodiment of Figure 15 can also advantageously be used to determine 
whether the ambient noise is primarily at 60 Hz, coiresponding to power line frequencies 
in the United States and Canada, or at 50 Hz, corresponding to power line frequencies in 
Europe. The foregoing modulation fi^uency of 316.7 Hz is selected to avoid the 
harmonics of the 60 Hz power line frequency as well as the 50 Hz power line frequency. 
If a significant shift in the power line frequency is detected such that aliasing of the 
ambient noise occurs at the firequencies of interest, then the modulation firquency can be 
changed to displace the modulation harmonics farther from the harmonics of the power 
fine frequency, such as, for example, by changing the 46,875 Hz sampling frequency, or 
by changing the modulus. 
Pre-Demodulation Decimation 

For convenience, the previous embodiments do not show the signal MF{K) being 
decimated before demodulation. However, as discussed in more detail below, the signal 
MFiK) can advantageously be decimated prior to demodulation. The pre-demodulation 
decimation technique can reduce the computational burden required to perform the 
donodulation operations, primarily because the decimated sample rate is lower than the 
original (undecimated) sample rate. Computation can also be reduced because, as will be 
seen, the numerical sequences used in the demodulator are, in some circumstances, shorter 
than the sequences given in Equations 32 and 33. Pre-demodulation decimation is a 
generalization of the previous embodiments and reduces to the previous embodunents 
when the pre-demodulation decimation rate is one. 

Figure 16 is a pictorial representation of a system that incorporates pre- 
demodulation filtering and decimation. Figure 16 is similar to Hgure 3, and like numbers 
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refer to like elements in the two figures. Figure 16 shows the first modulator 191 having a 
signal input SiO) and a modulation inpnt Afi(t). The second modulator 193 has a signal 
input Slit) and a modulation input ^2(0- The pair of signals S\{t) and S^t) represent the 
effect of the time-varying volume and scattering components of the blood in a fin^ (or 
other body part) on the red light and the inftared light, respectively, passing through tiie 
finger. The red light signal portion is caused by the variable attenuation of the red 
light passing through the finger 102 (shown in Figure 1). The inftared light signal portion 
52(0 is caused by tiie variable attenuation of the infrared light passing through the finger 
102. The ouQjuts of the first and second jiiodulators 191, 193 are provided to the receiving 
photodetector 150. The photodetector 150 is modded as an addw 194 and an adder 196. 
The outputs of the first and second modulators 191, 193 are provided to the addw 194 to 
genaate a composite signal M(t) where: 

M(0 = Si(/)A/i(0 + S2(t)MM (34) 
The ou^ut signal M(t) firom the adder 194 is provided to an adder 196 where a 
signal n(t) is added to the signal M(t). The signal n{f) represents a composite noise signal 
caused by ambient light (including DC and harmonics of &e power line frequency), 
electromagnetic picta^, and the like, which are also detected by the photodetector 150. In 
addition, the signal itQ) may also include noise at higher frequencies caused, for example, 
by other devices such as electrocauterization equipment, or die like. The ou^ut of the 
adder 196 is a signal Af{.t) = Mit) + n(t) which includes noise components as well as the 
signal components. 

The Af(0 signal ou^ut of the adder 196 (i.e., the output of the detector 150) is 
applied to the input of a signal processing bk>ck 1600. Within the signal processing blodc 
1600, the signal Af(t) is first passed through the amplifier 197 and then througji the analog 
ban^ass filter 198. The analog bandpass filter 198 provides anti-aliaang and removal of 
low frequency noise and DC. The filter 198 has a passband selected to pass dgnals in flie 
preferred range of 20 Hz to 10,000 Hz. The analog bandpass filter 198 removes a 
significant portion of the noise below 20 Hz. The signal components responsive to ttie 
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blood oxygen saturation are frequency shifted by the operation of the two modulation 
signals Miit) and and are passed by the analog bandpass filter 198. 

In one embodiment, the output of the analog bandpass filter 1 98 is sampled by the 
analog-to-digital converter 199 and converted therein to digital signals. In one 
embodiment, tiie signals are sampled at 46,875 samples per second. The digital signals 
from the analog-to-digital converter 199 are provided as inputs to a lowpass digital filter 
1620. Output signals bom the distal filter 1620 are provided to a sample rate 
compression block 1622 that reduces (compresses) tiie sample rate by a decimation rate Ry. 
The lowpass digital filter 1620 and sample rate compressor 1622 together comprise a 
decimator 1621 (decimation comprises lowpass filtering followed by sample rate 
compression). The digital filter 1620 provides anti-aliasing filtering and the sample rate 
compression block 1622 preferably operates at a sampling rate of at least twice tiie highest 
frequency of interest as determined by tiie digital filter 1620. In one embodiment, tiie 
sample rate compression block 1622 reduces tiie sample rate by a factor of R\ = 37, 
corresponding to tiie number of samples during tiie period x as illustrated in Figure 10. 
The output of the sample rate compression block 1622 provides one sample per time 
period T and tiius four samples per time period T. The output of tiie sample rate 
compression block 1622 is a signal MFik) (where A: is a discrete index) which comprises 
approximately 1,266 samples per second. 

The signal MFik) is provided as a first input to a first mixer 1624. The signal 
MF{k) is also provided as a first input to a second mixCT 1626. A first demodulating signal 
Di{k) is provided as a second iiq)ut to tiie first mixer 1624, and a second demodulating 
signal AC*) is provided as a second input to tiie second mixer 1626. The ou^ut of the first 
mixer 1624 is provided as an input to a first lowpass filter 1630, and tiie output of tiie 
second mixer 'is provided as an input to a second lowpass filter 1640. The bandwidths of 
tiie lowpass filteis 1630, 1640 are preferably approximately 10 Hz. The signal MF(k) is 
also provided as a firet input to a noise channel mixer 1628. A noise demodulating signal 
Do(ifc) is provided as a second input to tiie noise channel mixer 1628. The output of tiie 
noise channel mixer 1628 is provided to an input of a low pass filter 1650. The output of 
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the low pass filter 1650 is provided to a sample rate compressioii block' 1652. The output 
of the sample rate compression block 1652 is an estimate of flie noise n(t). The output of 
the lowpass filters 1630 is provided to an mpnt of a sample rate compressor 1632 and the 
output of the lowpass filter 1640 is provided to an input of a sample rate compressor 1642. 
5 The lowpass filter 1630 and the sample rate compressor 1632 together comprise a 
decimator 1631. The lowpass filter 1640 and the sample rate compressor 1642 together 
comprise a decimator 1641, 

The output of the decimator 1631 is a signal 5, (k) , which, as discussed below, is 

an estimate of the signal S\(k). The output of the decimator 1641 is a signal §2 (k) , which, 
10 as discussed below, is an estimate of the signal Siik). As will be shown below, the 
selection of the first demodulating signal D\{k) and the second demodulating signal DiQc) 
in accordance with the present invention can reduce or eliminate the effects of noise in the 

two output signals S^{k) and >$2(^) reduce or eliminate crosstalk between the 

two signals. 

IS The decimators 1632, 1642 decimate by a decimation rate Ri. In a prefeired 

CTibodiment, the decimators 1632, 1642 decimate by a decimation rate R2 = 20 to a 
sample rate of, for example, 633 Hz to provide a decimated ou^ut wldch can be fiirther 
processed in accordance with the methods and apparatuses described in the above- 
referenced patents. The decimations which occur in the decimators 1632, 1642 reduce the 

20 rate at which the output signals 5, {k) and S^ik) need to be processed v^le maintaining 
the sample rate well above the 10 Hz fi:equency content of the signals of interest The 
ou^uts of the decimators 1632, 1642 are provided on respective ou^ut lines 1634 and 
1644. 

Decimating the signal MF(k) prior to demodulation, altiiough not an approximation 
25 technique, can be simplified by assuming that each desired signal 5i(0 -does not change 
appreciably during each period t. In many applications it is reasonable to assume that the 
desired signals iS^i(0 and 52(0 ^ot change significanUy during the time int^val x shown 
in Figure 2. ^ One skilled in the art will recognize that a sufGcioit condition for this 
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assumption is that the highest significant frequency components in S\it) and Siit) are mudi 
lower than the modulation frequency. In the pulse-oximetry application the hig^ 
frequency of interest is typically around 10 Hz, which is far below the 316.7 Hz. 
fundamental of the modulation. Since n(/) is not a desired signal no such assumption is 
necessary for nit). Thus, while n{t) may vary erratically over a modulation cycle, the 
signals S\it) and 52(0 do not. Therefore, it is possible to perform pre-demodulation 
decimation that has UtUe effect on and 52(0 but may shape n(0 into n'(0. The 
measured signal is decimated by a factor /?i = Q (where g is the number of samples in a 
time period x) and then demodulated. 

Assuming R\ = Q, then the spectral domain representation of the signal A{F(k) at 
the output of the sample rate compression block 1622 is given by (approximately): 

w;4E[^,(/-7)v-7/s,(/.f)]^|t'.{/-f) (35) 

Since the sample rate compression block 1622 decimates at the same rate as the 
number of samples per period t, the decimation removes any x dependence in the 
expression for A£F(/). The frequency componaits indexed by m increase four times faster 
than the frequency components indexed by n. This occurs because the modulated signals 
Slit) and 52(0. which are indexed by n, occur in only one fourth of the samples, but the 
noise n{t), \Adch is indexed by m, occurs in eveiy sample. 

The demodulation operation can be perfbimed dther in the frequency or the time 
domain. A method for frequency domain demodulation of flie agnal A£F)[fc) can be 
obtained by rewriting Equation 35 as: 

MFO) A/F.,09+ MF_,0+A{F^0+ MF,(f)+ MF^(0+... (36) 
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where 

MF,W) = [5, (/) + "^i (/) + / T 

MF,(/) = [5,(/)-5,a)]/7' (37) 



Where n*(Jt) is the decimated noise signal n(t). Estimates for the signal can be 

obtained by shifting the spectra of MFi(/ ) and MF2(/) by -l/T and -2/7, respectively, 
and then dividing the sam of the resultant by 2. Likewise, ^iC/) can be obtained by 
dividing the differraice of the resultant spectra by 2. In other words: 

5,(/)=AiF,(/-l/D + A/Fj(/-2/r) 

Demodulation in the time domain is a more elegant method for obtaining Si(k) and 
Szik). Time domain demodulation is obtained by using the frequoicy shift property of the 
Foutier transfomi given by: 

F(c()+o>,) o e^**'/(0 (35) 

According to Equation 39, the frequency domain tenns ^lif) are related by a time shift 
in the time domiain and this property can be used to generate the demodulation sequences 
DtrDz. A more coiiq)lete development of this process (forthe general case of 7/channds) 
is provided in Equations 42-50 below and in the text acconqianying those equations. For 
the prcsHit case, where N-2, using equations 42-50 gives: 
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D^(k) = 0. 1, 0. 1. ... 
D,(k) = I. -0.5. 0. -0.5, ... (40) 
D^(k) = 0. -0.5. 1. -0.5, ... 

The sequences shown in Equation 40 are repeating sequences of the four values 
shown. Thus, the demodulation waveforms are no more than short repeating sequences of 
simple coefficients. Since the samples MF{k) are time domain sequences, demodulation 
simply involves multiplying the samples MF{]c) by the sequences in Equation 40. For 
example, the sequence of coefficients Dm = (0, 1, 0, 1, ...) is provided to the multipUer 
1628 to demodulate the signal AfFXJt) and produce the estimate of n(A). Similarly, the 
sequence of coefficients A(it) = (I, -0.5, 0, -0.5, ...) is provided to the multipUer 1624 to 
danodulate the signal MF{)C) and produce the estimate of 5i(*). 
Multiple Channel Modulation and D emodulation 

The two-channel pre-demodulation decimation technique described in the previous 
section can be extended to multi-channel systenas having more than two desired signals. 
Figure 17 illustrates an expansion of the two-channel modulator into a multi-channd 
modulator/demodulator. Figure 17 shows the first modulator 191 and the second 
modulator 193 as shown in Figure 16. Further, Figure 17 shows a third moduhitor 1701 
and an A^* modulator 1702. The signal mput S\{f) and a modulation input M\{i) are 
provided to the first modulator 191. The signal input Sz{() and a modulation input Mrii) 
are provided to the second modulator 193. A signal input 53(0 and a modulation input 
A/3(0 are provided to the third modulator 1701. A signal input S«(0 and a modulatioa 
inputM<0 are provided to the modulator 1702. 

The photodetector 150 is modeled as an adder 194 and an adder 196. The outputs 
of the modulators 191, 193, 1701. and 1703 are added together in the adder 194, to 
genoate a composite signal A/(0 where: 

M(0 = 5,(0A^.(0 + 5,(/)A/,(0 + 5,(0A/,(0+...+5^(0A/|,(0 (41) 
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The signal M(t) from the adder 194 is provided to the adder 196 where the signal 
M(t) is added to the signal nit) which rq>resaits a compoate noise signal caused by 
ambient light, electromagnetic pickup, and the like, y/bich are also detected by the 
photodetector 150. The output of the adder 196 is the signal M(t) = A/(/) + n(0, v^uch 
includes the noise components as well as the signal components. 

The Afit) signal output of the adder 196 (i.e., the output of the detector 150) is 
applied to the input of the signal-processing block 1700. Wthin the signal-processing 
block 1700, the signal A/(0 is first passed through an amplifier 197 and then through the 
analog banc^ass filter 198. The analog bandpass filto: 198 provides anti-aliasing and 
ronoval of low fiequoicy noise and DC. The desired signal components in the signals 
S^t) are frequoicy ^fted by the operation of the modulation signals M/(t) and are passed 
by the analog bandpass filter 198. 

The output of the analog bandpass filter 198 is sanpled by the analog-to-digital 
converter 199 and converted therein to digital signals and provided to an input of the 
lowpass digital filter 1620. Ou^ signals fiwm the digital filter 1620 are provided to a 
sample rate compression block 1622, which reduces the sample rate by a decimation factor 
^1. Together, the distal filter 1620 and the sample rate compression block 1622 comprise 
a decimator 1621. The output of the sample rate compression block 1622 is a signal 
MI\k). The signal MFlk) is provided as: the first input to the first mixer 1624; the first 
input to the second mixer 1626; a first input to a third mixer 1710; a first input to an 1th 
mixer 1712; and a first input to a noise channel mixer 1713. A first demodulating signal 
Diik) is provided as a second iiq>ut to the first mixo^ 1624. A second demodulating agnal 
Diik) is provided as a second input to the second mixer 1626. A third demodulating signal 
Diik) is provided to the third mixer 1710. A fourth demodulating signal D^k) is provided 
to the iV* mixer 1712. A noise demodulating signal Ddk) is provided to the noise channel 
mixer 1713. The outputs of the mixers 1624, 1626, 1710. 1712, and 1713 are provided as 
respective inputs of the lowpass filters 1630, 1640, 1720, 1730, and 1740. The outputs of 
the lowpass filters 1630, 1640, 1720, 1730, and 1740 are provided as respective inputs of 
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the decimatois 1632, 1642, 1721, 1731 and 1741. Each of the declinators 1632, 1642, 
1721, 1731 and 1741 reduces the sample rate by a decimation rate Ri. 

The ou^ut of the sample rate compressor 1632 is a signal 5,(Ar), which, as 
discussed below, is an estimate of the signal S\{k). Likewise, the output of the sample rate 
compressor 1642 is an estimate of 52(0. the output of the sanqjle rate compressor 1721 is 
an estimate of the signal S^it), the output of the sample rate compressor 173 1 is an estimate 
of the signal 5n(/), and the ou^ut of the sample rate compressor 1741 is an estimate of the 
signal n(0- 

As will be shown below, the selection of the donodulating signals D^t) for i - 
Q..I/ in accordance with die present invention can substantially reduce or eliminate the 
effects of noise in the output signals Si(k) and and can also substantially reduce or 
eliminate crosstalk between the signals. 

As shown in Figure 17, a set ofN+l signals Si[k] and n(k) are sampled at 

a rate T7QN, where 7 is a modulation period For simplicity, the decimation rate R\ is 
assumed to be the same as the factor Q. The assumption that i?i = Q is not a necessary 
assumption, but rather is used h«e to sin5)lify the mathematics. The signals are combined 
according to the formula: 

S(k) = M,(k)S,(k)+ M2(k)Si(k)+ Mj(k)Ss(k)+ ... + MN(k)SN(k)+n(k) (42) 
Using the symbol * to denote the convolution operator, the terms M,{k) are given by: 

A/,(A) = A(2M/r)*P,(0 Lr/oA, 

M,(*) = A(2M/r)*P,(0L„c,. ^^^^ 
M,{k) = A(2M IT)*P^{t) 



where 



[0 othennse 
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and 



Pi(O^Y. ^(^'""^ (45) 

(where b{k) is the Kroneker delta function, which is 1 for /: = 0, and 0 for all other values 
of A), and 



Pi(y = Pi(/-~^) for i^2...N (46) 



After the pre-demodulation and sample rate compression stage 1622, which decimates by a 
factor the signal in the frequency domain is given approximately by 



AfF(/)4£[5.(/-^).4:5,(/-f)^.^;..^,(/-f)] 



(47) 



where 



4; = (48) 



The demodulator sequences are then given by: 
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(49) 

where 

P,(k) = P,(t)U}L. P:(k) = P2(0U^ PH(k) = P„(t)\.^ (50) 

IN 2« 

The post demodulation lowpass filters 1630, 1640, 1720, 1730 and 1740, and the 
post demodulation sample rate compression stages 1632, 1642, 1721, 1731 and 1741 
suppress high frequency artifacts which are produced by the modulation/demodulation 
process. Note that Equation 49 reduces to Equation 40 for //= 2. 
Adaptive Demodulation 

The multi-channel pre-demodulation decimation technique described in the 
previous section can be extended to an adaptive multi-channel system having an adjustable 
pre-demodulation decimation rate and an adjustable post-demodulation decimation rate. 
Figure 18 illustrates an expansiim of tiie multi-channel modulator into a adaptive multi- 
channel modulatoi/demodulator 1800. Figure 18 shows the first modulator 191 and the N 
* modulator 1702 as shown in Figure 17. The signal input Si(t) and a modulation input 

are provided to the first modulator 191. A signal input Svit) and a modulation input 
Mail) are provided to the //* modulator 1 702. 
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The photodetector 150 is modeled as an adder 194 and an add^ 196. The outputs 
of the modulators 191, 193, 1701, and 1703 are added together in the adder 194, to 
generate a composite signal M{t) where: 

M(0 = 5,(OM,(0+...+SK(OAfw(0 (51) 

The signal M(0 &om the adder 194 is provided to the adder 196 where the signal 
M(t) is added to the signal n(0 which represents a composite noise signal caused by 
ambient light, electromagnetic pickup, and the like, \»*ich are also detected by the 
photodetector 150. The output of the adder 196 is the signal (0 = Mit) + nit), wMdi 
includes noise components as well as the signal components. 

The AfQ) signal output of the adder 196 (uc^ the output of the detector 150) is 
applied to the input of the signal-processing block 1800. Wthin the dgnal-processing 
block 1800, the signal (r) is first passed through the amplifier 197 and then through the 
analog bandpass filter 198. The analog bandpass filter 198 provides anti-aUasing and 
removal of low frequency noise and DC. The desired signal components in the signals 
are frequency shifted by the operation of fbs modulation signals M^t) and are passed 
by the analog bandpass filter 198. 

The ou^ut of the analog bandpass filt» 198 is sampled by the analog-to-digital 
converter 199 and converted therein to digital signals and provided to an input of a 
decimation block 1820. The adaptive decimation block 1820 comprises a digital lowpass 
filt» and a sanq)le rate compressor that reduces Hbs sample rate by the decimation rate Ri. 
The filter coefficients and decimation rate i?i are provided to a control mput of the adap6vc 
decimation block 1820 by an output of an adaptive algorithm block 1850. Equation 35 
assumes that the decimation rate /?, is equal to Q. However, in general, the vatae of Q 
may be different than the decimation rate iZj. The output of the adaptive decimation blodc 
1820 is a signal AfFX*). 

The signal MFlk) is provided to the first input of the first mixer 1624, to the first 
it^ut of the th mixer 1712, and to the first input of the noise channel mixer 1713. A 
first demodulating signal Di(k) is provided to a second input of ttie first mixer 1624 from a 
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signal generator 1841. The fourth demodulating signal D^k) is provided to the //* mixer 
1712 from an output of a signal generator 1831. The noise demodulating signal Do(A) is 
provided to the noise channel mixer 1713 from an output of a signal generator 1832. A 
control input to each of the signal generators 1831, 1832, and 1841 is provided by the 
output of the adaptive algorithm 1850. In yet another embodiment, the adaptive algorithm 
1850 may also be controlled by other signal processing elements downstream of the signal 
processor 1800. 

The outputs of the mixeR 1713, 1624, and 1712 are provided as respective inputs 
to adaptive decimation blocks 1840, 1830, and 1834 respectively. Each of the adaptive 
decimation blocks 1840, 1830, and 1834 has a control input provided by the ou^ut of the 
adaptive algorithm block 1850. The output of the adaptive decimation block 1840 is an 
estimate of the signal n(t) and it is provided to an input of the adaptive algorithm block 
1850. In an alternate embodiment, the signal estimates S,ik) are also provided to the 
adaptive algorithm block 1850. 

An output of the decimator 1830 is a signal , which, as discussed above, is 
an estimate of the signal S\(k). likewise, the output of the decimation block 1834 is an 
estimate of the signal S^t). As shown above, the sdection of the demodulating signals 
A<0 for i = O..Jf in accordance witfi the present invention substantially reduces or 
eliminates the effects of noise in the output signals and n{kX and also substantially 
reduces or eliminates crosstalk between the signals. 

As shovm in Figure 18, a set o(N+l signals Si[k] i = 1...//, and nik) are sampled at 
a rate I7QN, where T is a modulation period, and R\ is the decimation rate of the 
dedmation block 1820. The signals are combined according to the formula: 

S(k)='Mi(k)S,(k)+...+MN(k)Ss(k)+n(k) (52) 

Each of the adaptive decimators 1820, 1840, 1830, and 1834 comprises a digital 
lowpass filter and a sample rate compressor. The characteristics of flie digital lowpass 
filters (e.g., the number of filter coefiicients and values of the filter coefficients) and the 
sample rate compression factor of each adaptive decimator is provided to a control input of 

-45- 



Substitute Sheet (Rule 26) 



o o 



wo 99/52420 PCT/US99/07825 



tte adaptive declinator. The control inputs are driven by an adaptive algorithm 1850. The 
signal generators 1831, 1832 and 1841 generate the demodulation sequences for the 
demodulators 1624, 1712, and 1713 respectively. The demodulation sequences produced 
by the signal generators 1831, 1832 and 1841 are controUed by the adaptive algorithm 
1850. 

The adaptive algorithm adjusts the pre-demodulation decimation rate i2i (in the 
adaptive demodulator 1820), and the post-demodulation decimation rate R2 (in the 
adaptive demodulators 1830, 1834 and 1840) accordmg to tiie noise in die noise estimate 
n{k) 1746 and (optionally) according to Uxe signals S,ik) . The product RiR2 is the total 
decimation rate fiom the signal S(k) at the output of the A/D converter 199 to the signals 
5.(ifc)at the output of the signal processing block 1800. The adaptive algorithm may 
adjust R\ and Ri such that the product R\R2 varies, or the adaptive algorithm may adjust Ri 
and R2 such that the product R\Ri is substantially constant. Typically, the adaptive 
algorithm will keep the RxRi product constant so that the signal procesang blodcs 
downstream of the signal processor 1800 will operate at a substantially constant sample 
rate. 

Typically, each of Ae signal generators 1841, 1831 and 1832 generates a repeating 
sequence of numbers. The number of elements in tiie sequaice is a fimction of ttie 
decimation fectoriii. As discussed above in connection with Figure 3, \i4ieni?i =1, there 
are preferably 148 values in eadi demodulation sequence. As discussed above in 
connection witii Figure 17, when Ri = 37, there arc preferably only 4 values in the 
demodulation sequences. 

The adaptive algorithm selects Ru Rz, and the filter transfer functions in the 
adaptive deciihators 1820, 1830, 1834, and 1840 to improve the quality of ±e output 
signals S,ik). For example, in hi^ ambient noise environments, tiie higher orda: 
hamionics of tiie output signals are often contaminated by ambient noise (as discussed in 
connectim wifli Figures 14 and 20). Thus, tiw higbCT order harmonics are prrferably not 
demodulated when ambient noise is present To avoid demod u lati o n of the higher order 
harmonics the adaptive demodulator 1850 can set /?i = 1 and ii2 = 37, and thereby 
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demodulate accoiding to the method described in comiection with Figures 3-14. 
Alternatively, the adaptive demodulator 1850 can set R\ = 37, set /?2 = U and set the 
transfer function of the lowpass filter in the adaptive decimator 1820 to provide a very fzst 
rolloff (thereby filtering out the higher order hannonics). 

Conversely, in low ambient noise envirorunents, the higher order harmonics of the 
output signal are less contaminated by ambient noise, and thus the higher order harmonics 
may be demodulated In one embodiment, to demodulate the higher order hannonics, the 
adaptive demodulator 1850 can set R\ = 37 and set iJj = U to demodulate according to the 
method described in connection with Figure 17. This is especially advantageous when 
perfusion is low, because, when perfiision is low the ou^ut signals S^ik) are typically 
very weak and are contaminated by random noise. Demodulating more of the higher order 
harmonics increases the signal-to-noise ratio because it adds the harmonics (which are 
correlated) to the output signals, and tends to average out the noise (which is uncorrelated). 
Thus, the signal strength increases, and the noise is reduced. 

One skilled in the art will recognize that the examples in the preceding two 
paragraphs are merely two points on a continuum and that the adaptive algorithm 1850 can 
generate many desirable solutions on the continuum. 
Ambient Light Rejection 

In the pulse oximeter, one of the major contributors to the noise signal n(t) is 
ambient light that is detected by the photodetector 150. One aspect of the present 
invention advantageously provides a method for choosing the modulation sampling lateyj 
and the factor Q so that the effects of ambient light can be removed by the post 
demodulation filtering and decimation stages. Note that g is the number of samples 
during the on.period C-e., modulation signal sample turn on time Q) and is preferably also 
the decimation rate R\ for the pre-demodulation sample rate compressor 1622 (in general 
the values of Q and R\ may be different). The particular embodiment described by 
Equation 35 assumes that the value Q is also used as decimation rate R\ for the pre- 
demodulation decimator 1820. 
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bi the system diown in ngmes 3 and 16, which demodulates two hamionics, the 
period of a modulation cycle is givoi by: 



T = 4Q/f, 

where^ is the sample rate. Defining the two line equations 



(53) 



where 



- line frequencies of concern 
n = line frequency hannonic numbers of concern 



(54) 



(55) 



then the effects due to ambient light will be nunimized when 

\y{f^. ny^SBF 
\z{f,.n)'^SBF 



(56) 



where SBF is the stop band frequency of the post demodulation and decimation stages 
(e.g., the 10 Hz iowpass filter 1630 and the sample rate compressor 1632, etc.). 

Rgure 19 is a flowchart showing a m^od for selecting^ and Q. The method 
begins at a process block 1902 viierein the ambient light firequencies^ and important 
hannonic conqranents n are identified. Inqwrtant harmonics are defined as those 
harm<»iics that will degrade system pofinmance below accqytable levels ^en detected by 
the detector 150. The process then advances from the process block 1902 to a process 
block 1904. hi the process blodc 1904, the values o^f^ and n identified in the process 
Mode 1902 are used in corQunctim with Equation 54 to identify a coUectitm of accqrtable 
values of T. Upon completion of the process block 1904, the process advances to a 
process block 1806. In the process 1906, suitable values of^ and Q are diosen uang tbe 
values of T obtained in the process block 1904 and the equation r= AQIf^. One skilled m 
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the art will recognize. that, since T is proportional to the ratio of Qlf^^ knowing Twill not 
uniquely determine either ^ or Q. 

For example, given power line frequencies of 50 ± 1 Hz and 60 ± I Hz then the 
range o{f^ is given by approximately the union of the interval 49-51 Hz and the interval 
59-61 Hz, which can be expressed mathematically as: 

^ [49.51] KJ [59.61] (57) 

Assuming that all harmonics up to the 18* harmonic are to be suppressed, then n = 1^.18. 
In a preferred embodiment, using these values for^ and n, application of the method in 
Figure 19 results in^s = 46,875 Hz and acceptable Q values of 37 and 41. 

The process leading to Equation 57 is illustrated graphically by Figure 20, vAiere 
the harmonics of the ambient light frequency (in Hz) are plotted versus the 
plethysmograph signal frequency (also in Hz). Figure 20 has an x-axis showing the 
ambient light frequency from 44 Hz to 64 Hz. The ambient light frequency will usually 
correspond to the frequency of the power lines, which is nominally 60 Hz (in the U.SO 
and 50 Hz (outside the U.S.). However, power line frequency regulation typically 
varies somewhat, and thus Figure 20 shows frequencies above and below the nominal 
frequencies. 

Figure 20 also shows a y-axis showing the plethysmograph signal frequency from 
-10 Hz to 10 Hz. One skilled in the art will recognize that negative frequencies occur in 
the mathematics described above. In particular, a signal that is modulated from baseband 
up to some carrier frequency will exhibit two sidebands, a sideband above the carrier 
frequency corresponding to the frequency of the baseband signal, and a adeband below the 
carrier ftequepcy corresponding to the negative of the baseband frequency. Thus, ^^Aum 
dealing with modulation and demodulation, it is convenient to deal with positive and 
negative frequencies. 

Figure 20 also shows harmonic lines corresponding to the 5*, 6*, 7*, 10^, 11*, 
12*, 13*, and 14* harmonics of tiie ambient light frequency. The harmonic lines 
correspond to tiie harmonics produced in the plethysmograph signal by the demodulation 
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(mixing down) of hannonics of fhe power line frequency. The lines in Figuxe 20 are 
calculated using Equation 54 for 1/7- 316.72 Hz. Some of the harmonic lines conespond 
to y{f^) and some correspond to z(/^) from Equation 54. EUomonic lines that are not 
shown (e.g., the line corresponding to the 8^ harmonic) M outside the displayed limits of 
the X-axis and y-axis. 

Figure 20 can be used to detennine the stop band frequencies as shown in Equation 
56. For example, the harmonic lines in Figure 20 show that for an ambient light frequency 
of 49 Hz, the 13* harmonic of the ambient light frequency will appear in the 
plethysmograph signal at approximately 3 Hz. Thus, Figure 20 shows that for 
plethysmograph bandwidth of 10 Hz, none of the first 14 harmonics of the ambient light 
will appear in the plethysmograph signal for ambient light frequencies between 
approximately 612 Hz and approximately 58.5 Hz» which is consistent with Equation 57. 
The first ambient harmonics that do appear for a plediysmogRq)h bandwidth of 10 Hz are 
the 5* harmonic and the 1 1* hannonic. 
Other Embodiments 

In the preferred embodunent of the present invention, the hardware described 
above is implemented in a digital signal processor and associated circuitry. The LED 
modulation block 104 and the LED dCTiodulation state table block 352 comprise 
algorithms implemented by program code executed by the digital signal processor. In 
addition, the configuration variables, such as for example, the hardware delay value, the 
iKUxiware distortion value and the hardware scale value are provided as inputs to the 
digital signal processor when it is set up. For example, the main operating program of 
the digital signal processor may be stored in non-volatile ROM or PROM, and the 
variables may be stored in flash memory during a setup procedure. Tedmiques for 
communicating to and from a digital signal processor during sudi setup procedures are 
well known to persons of skill in the art, and will not be desoibed in detail heroin. For 
example, the configuration bus 310, discussed above, represents a communication path 
to the flash mmory during such a setup procedure. The data provided to the 
configuration bus 310 may be provided by a system operator (not shown) or the data 
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may be provided from look-up tables (not shown) maintained for different embodiments 
of the LEDs 106, 108 and the detector ISO. 

Although described above in connection with a pulse oximetry system wherein a 
parameter to be measured is the attenuation of red and infiared light passing through a 
portion of a subject's body, it should be understood that the method and apparatus 
described hsrein can also be used for other measurements where two or more signals are 
passed through a system to be analyzed. In particular, the present invention can be used to 
demodulate two combined parametric signals responsive to the system to be analyzed 
where the two parametric signals have a predetermined timing relationship between them, 
as described herein. 

One skilled in the art will recognize that the lowpass filters provided in connection 
with the decimation blocks may provide other filter fimctions in addition to lowpass 
filtering. Thus, for example, the lowpass filters 1620, 1622, 1630, 1640, 1650, 1720, 
1730, and 1740, and the decimators 1820, 1830, 1834, and 1840 may provide other filter 
functions (in addition to lowpass filtering) such as, for acample, bandpass filtering, 
bandstop filtering, etc. Moreover, the post-demodulation decimation rate Ri need not be 
the same for each output channel. Thus, for example, in Figure 18, the decimator 1840 
may have a first decimation rate Ri = n while the decimators 1830 and 1834 have a second 

decimation rate /l2= ^2. 

Although described above in connection with a particular embodunent of the 
present invention, it should be understood die description of tfxe embodiment is illustrative 
of the invention and are not intaidcd to be limiting. Various modifications and 
plications may occur to those skilled in the art without departing fiwm die tme spirit and 
scope of the invention as defined in the s^pended claims. 



-51- 



Substitute Sheet (Rule 26) 



o o 

wo 99/52420 PCT/US99/07825 

WHAT TS CLAIMED IS : 

1 1 . An apparatus for measuring blood oxygenation in a subject, said apparatus 

2 comprising: 

3 a first signal source which zpp&es a first input signal during a first time 

4 interval; 

5 a second signal source which applies a second input signal during a second 

6 time interval; 

7 a detector which detects a first parametric signal responsive to said first 

8 input signal passing through a portion of said subject having blood therein and 

9 whidi detects a second parametric signal responsive to said second input signal 

10 passing through said portion of said subject, said detector generating a detector 

11 output signal responsive to said first and second parametric signals; and 

12 a signal processor which receives said detector ou^ut signal, said signal 

13 processor demodulating said detector output signal by applying a first 

14 demodulation signal to a signal responsive to said detector ou^)ut signal to 

15 generate a first output signal responsive to said first parametric signal and applying 

16 a second demodulation signal to said signal responsive to said detector output 

17 signal to generate a second output signal responsive to said second parametric 

18 signal, each of said first demodulation signal and said second demodulation signal 

19 comprising at least a first component having a first firequency and a first amplitude 

20 and a second component having a second frequency and a second anq}litude, said 

21 second frequoicy being a harmonic of said first firequency, said second amplitude 

22 selected to be related to said first anoqplitude to minimize crosstalk firom said first 

23 paramebic signal to said second ou^t ^gnal and to minimize ax>ssta]k fix>m said 

24 second parametric signal to said first output signal 

1 2. The mettiod of Claim 1, ^odn the second amplitude is detemtiined by 

2 turning off one of the first and second signal sources and measuring the crosstalk between 

3 one of the parametric signals and the non-corresponding output signal vMe varying the 

4 second amplitude and selecting a second amplitude which minimizes tiie measured 

5 crosstalk. 
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3. A method of miimnizing crosstalk between two stgoals generated by 
applying a first pulse and a second pulse to measure a parameter, wherein said first pulse 
and said second pulse are applied periodically at a first repetition rate defining a period, 
and wherein said first pulse is generated during a first interval in each period and said 
second pulse is generated during a second uiterval in each period, said second interval 
spaced apart ftom said first interval, said first and second pulses producing first and second 
parametric signals responsive to said parameter, said first and second parametric signals 
being received by a single detector which outputs a composite signal responsive to said 
first and second parametric signals, said ntethod comprising the steps of: 

applying a first demodulation signal to said composite signal to generate a 
first demodulated ou^ut signal, said first demodulation signal comprising at least a 
first component having a first fi:equcncy corresponding to said first repetition rate 
and having a first amplitude, said first demodulation signal fiirther comprising a 
second component having a second firequency which is a harmonic of said first 
fiequency and having a second amplitude which has a selected proportional 
relationship to said first amplitude; 

applying a second demodulation signal to said composite signal to generate 
a second demodulated ou^ut signal, said second demodulation signal comprising 
said first component at said first fi:cquency and said first amplitude and comprising 
said second component at said second firequency and said second amplitude, at 
least one of said first and second components of said second danodulation signal 
having a selected phase difference with respect to the corresponding one of said 
first and second components of said first demodulation signal; and 

lowpass filtering said first demodulated output signal to generate a first 
recovered ou^ut signal responsive to said first parametric signal; and 

lowpass filtering said second demodulated output signal to generate a 
second recovered ou^ut signal responsive to said second parametric signal 
4. The method as defined in Claim 3, wherein said selected phase difference 

iSTC. 
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5. The method as defined in Claim 3, wherein: 

said first pulse and said second pulse aie geneially rectangular pulses 
having a duty cycle, and wherein said rectangular pulses comprise a plurality of 
sinusoidal components includii^ a fimdamental conqx>nent corresponding to said 
first frequency and a first harmonic conqx)nent corresponding to said second 
frequency, said fimdamental component having a fimdammtal component 
amplitude and said first harmonic conq)onent having a first harmonic component 
amplitude, said first harmonic component amplitude related to said harmonic 
component amplitude by a first proportionality value; and 

said second amplitude of said second component of said first demodxdation 
signal is related to said first amplitude of said first component of said first 
demodulation signal by a second proportionality value which is approximately the 
inverse of said first proportionality value. 

6. The method as defined in Claim 3, fiirther includmg the steps of: 
sampling said composite signal when neither said first pulse nor said 

second pulse is active to obtam a sampled signal; and 

measuring said sampled signal to determine a noise level of said parametric 

signals. 

7. The method as defined in Claim 3, fiirther including the steps of: 
perfomiing a transform on said composite signal to generate a spectra of 

said composite signal; 

sampling said spectra at a plurality of frequencies other than at 
predetermined ranges of frequencies around said first frequency and around 
harmonics of said first firequeacy; 

determining an average of the magnitudes of said sampled plurality of 
fi:equencies; and 

comparing said average to a selected threshold to determine whether the 
avemge m^nitude »ceeds said selected tiireshold 
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g. A method of demodulating a composite signal generated by applying fixst 
and second periodic pulses of electromagnetic energy to a system having a parameter to be 
measured and by receiving signals responsive to said electromagnetic energy afker having 
passed through said system and being affected by said parameter being measured, said 
signals received as a composite signal having components responsive to said first and 
second pulses, said method comprising the steps of: 

applying a first demodulation signal to said composite signal to generate a 
first demodulated signal, said first demodulation signal comprising a first 
component having a first fiwjuenCy corresponding to a repetition frequency of said 
first and second pulses and comprising a second component having a fi^uency 
which is a harmonic of said first firequency, said first component having a first 
amplitude and said second component having a second amplitude, said second 
amplitude having a predetermined relationship to said first amplitude, said 
predetemuned relationship selected to cause said first donodulated signal to have 
low firequency components responsive only to said first pulse; and 

lowpass filtering said first demodulated signal to generate a first output 
signal, said first output signal varying in response to an effect of said parameter on 
the electromagnetic energy received fix3m said first pulse. 
9. The method as defined in Claim 8, further including the steps of: 

applying a second demodulation signal to said composite signal to generate 
a second demodulated signal, said second demodulation signal having first and 
second components corre^nding to said first and second components of said first 
demoduladon signal, at least one of said first and second components of said 
second demodulation signal having a selected phase relationship with the 
corresponding one of said first and second components of said first demodulation 
signal; and 

lowpass filtoing said second demodulated signal to generate a second 
output signal, said second ou^ut signal varying in response to an effect of said 
parameter on the electromagnetic energy received fixmi said second pulse. 
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1 10. The method as defined in Claim 9, wherein said select^ 

2 isaTcphasedifiference. 

1 1 1 . A pulse oximetiy system, comprising: 

2 a modulation signal generator, said modulation signal generator generating 

3 a first modulation signal comprising a first pulse which repeats at a first repetition 

4 frequency, said first pulse having a duty cycle of less than 50%, said modulation 

5 signal generator generating a second modulation signal comprising a second pulse 

6 which also repeats at said first repetition frequency, said second pulse having a 

7 duty cycle of less than 50%, said'second pulse occurring at non-overls^ping times 
6 with respect to said first pulse, said first and second pulses comprising a plurality 
9 of components wherein a first component has a frequency corresponding to said 

10 repetition firequency and a second c(Hnponent has a second fi:cquency 

11 corresponding to twice said first frequency, said second component having an 

12 amplitude which has a first predetermined relationship to an amplitude of said first 

13 component; 

14 a first transmitter which emits electromagnetic energy at a first wavelength 

15 in response to said first pulse; 

16 a second transmitter which emits electromagnetic energy at a second 

17 wavelength in response to said second pulse; 

18 a detector which receives electromagnetic energy at said first and second 

19 wavelengths after passing through a portion of a subject and whidi genoiates a 

20 detector output signal responsive to the received electromagnetic energy, said 

21 detector ou^ut signal includii^ a ^gnal component responsive to attenuation of 

22 said electromagnetic energy at said first wavelength and a signal conqxmcot 

23 responsive to attenuation of said electrom^etic energy at said second 

24 wavelengdi; 

25 a first dCTiodulator which multiplies said detector signal by a first 

26 demodulation signal and generates a first demodulated output signal, said first 

27 demodulation signal comprising a first component having said first frequency and 
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28 having a first anq)Iitude and comprising a second component having said second 

29 frequency and having a second amplitude, said second amplitude having a second 

30 predetemuned relationship to said first amplitude vMch second pred^ennined 

31 relationship is inversely proportional to said first ptedetenmned relationship; and 

32 a second demodulator viiich multiplies said detector signal by a second 

33 demodulation signal and generates a second demodulated output signal, said 

34 second demodulation signal comprising a first component havhig said first 

35 frequency and having said first amplitude and comprising a second component 

36 having said second fiequency and having said second amplitude, at least one 

37 component of said second demodulation signal having a selected phase 

38 relationship with a correspondii^ one component of said first demodulation signaL 

1 12. The method as defined in Claim 11, wh^ein said selected phase 

2 relationship is a tc phase difference* 

1 13. A method of minimizing crosstalk between two signals generated by 



2 applying a first pulse and a second pulse to measure a parameter, wherein said first pulse 

3 and said second pulse are applied periodically at a fiirst repetition rate defining a period, 

4 and wh^ein said first pulse is generated during a first interval in each period and said 

5 second pulse is generated during a second interval in each period, said second interval 

6 spaced apart torn said first interval, said first and second pulses producing first and second 

7 parametric signals responsive to said parameter, said first and second parametric signals 

8 being received by a single detector which outputs a composite signal responsive to said 

9 first and second parametric signals, said method comprising the steps of: 



10 providing said composite signal to an analog to digital converter to produce 

11 a sequence of digital values; 

12 decimating said sequoice of digital values to produce a decimated 

13 sequence of digital values; 

14 applying a first sequence of demodulation coefficients to said decimated 

15 sequence of digital values to gmerate a first donodulated output signal; 
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applying a second sequence of demodulation coefficients to said decimated 
sequence of digital values to generate a second demodulated ou^iut signal; 

lowpass filtering said first demodulated output signal to generate a first 
recovered ou^ut signal responsive to said first param^c signal; and 

iowpass filtering said second demodulated ou^ut signal to generate a 
second recovered output signal responsive to said second parametric signal. 

14. The method as defined in Claim 13, further comprising the stq)s of: 
applying a third sequence of demodulation coefficients to said decimated 

sequence of digital values to generate a third demodulated output signal; and 

Iowpass filtering said third donodulated output signal to generate a 
recovered output signal responsive to noise detected by said detector. 

15. A method of demodulating a composite signal gen^ed by applying a 
plurality of periodic pulse trains of electromagnetic energy to a system having a parameter 
to be measured and by receiving signals responsive to said electromagnetic energy after 
having passed through said system and being affected by said parameter being measured, 
said signals received as a composite signal having componoits responsive to said plurality 
of periodic pulse trains, each of said components corresponding to one each of said pulse 
trains, said method comprising the steps of: 

sampling said composite signal using an analog to digital converter to 
produce a sequence of digital values; 

decimating said sequence of digital values to pnxiuce a decimated 
sequence of values; 

providing said decimated sequmce of values to a first ii^ut of each of 
phirality of multipliers; 

providing a plurality of sequences of danoduladon coeffidents to a second 
iiq>ut of each of said multipliets sudi that each multiplier is provided a unique 
sequence of demodulation coeffidrats; and 

Iowpass filtering the output of each of said multipliers such that the output 
signal of each of said Iowpass Slters approximately corresponds to one each of said 
components. 
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1 16. A pulse oximetiy system, comprising: 

2 a modulation signal generator, said modulation signal generator generating 



a first modulation signal comprising a first pulse which repeats at a first repetition 
fircquency, said first pulse having a duty cycle of less than 50%, said modulation 
signal generator generating a second modulation signal comprising a second pulse 
which also repeats at said first repetition frequency, said second pulse having said 
duty cycle of less than 50%, said second pulse occurring at non-overlapping times 

8 with respect to said first pulse; 

9 a first transmitter vduchcrnitselectrornagnetic energy at a first wavelength 
io in response to said first pulse; 

^1 a second transmitter vdiich emits electromagnetic energy at a second 

^2 wavelength in response to said second pulse; 

^3 a detector which receives electromagnetic energy at said first and second 

wavelengths after passing through a portion of a subject and which generates a 
detector output signal responsive to the received electromagnetic en^gy, said 
16 detector output signal including a signal componmt responsive to attenuation of 

said electromagnetic energy at said first wavelength and a signal componait 
responsive to attenuation of said electromagnetic energy at said second 

19 wavelength; 

20 a sampling analog to digital converter which converts said detector output 

21 into a sequence of digital values; 

22 adecimator which decimates said sequence of digital values to produce a 

23 decunated sequence; 

24 a first demodulator which multiplies said decimated sequence by a 

25 demodulation sequence and generates a first demodulated output signal; and 
a second demodulator which multiplies said decimated sequence by a 

second demodulation sequence and generates a second demodulated output signal 
17. . The apparatus of claim 16 fiirther comprising a third demodulator which 
2 multipUes said decimated sequence by a third demodulation sequence and generates a 
ttiird demodulated output signal corresponding to noise produced by said detector. 
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1 18. The iq)paratus of claim 16 wherein said decimator hasr a decimation rate 

2 equal to the number of said samples produced during said duty cycle* 

1 19. In a systffli that includes: a sampling fipequency generator generating a 

2 sampling frequency^; a detector; a modulation signal generator, said modulation signal 

3 generator generating a sequence of pulses having a puke repetition frequency, said pulses 

4 having a duty cycle of Q sample periods of said sampling frequency yi; a transmitter which 

5 emits electromagnetic energy at a wavelength in response to said pulse, said detector 

6 receiving said electromagnetic energy to generate a detector output signal responsive to the 

7 received electromagnetic energy, said detector output signal including noise caused by 

8 ambient electromagnetic energy detected by said detector; a digital to analog convertar 

9 which produces distal sanq)les of said detector output signal at said sampling frequency; 

10 and a demodulator which demodulates said digital samples to produce a desired output 



11 signal, a method for m^^t^r^^^^g noise due to ambient electroms^etic energy in said 

12 desired output signal comprising the steps of: 

13 identifying all undesired frequency con^xments of said ambient 

14 electromagnetic energy detected by said detector, 

15 using said imdesiied frequmcy components to compute a set of acceptable 

16 modulation cycle times r, and 

17 using said acceptable modulation cycle times to select saidy^ and said Q 

18 using an equation T=4Q/fs. 

1 20. A method for demodulating a multi-channel composite signal comprising 

2 the acts of: 

3 generating a demodulation sequence for a selected chaimel; 

4 providing said demodulation sequence to a first input of a demodulator, 

5 providing a sanq^led conq^osite signal to a second ixsput of said 

6 demodulator 

1 21. The method of Claim 20, further comprising the act of decimating an 

2 ou^ut of said demodulator by a decimation factor R. 

1 22. Themethodof Oaim 21, wherein said act of deciniating includes the acts 

2 of: 
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lowpass filtering; and 
sample rate compressing. 

23. The method of Clahn 20, further comprising the act of decimating said 
sampled composite signal before providing said sampled composite signal to said second 
input of said demodulator. 

24. The method of Claim 23, wherein said act of decimating includes the acts 
of lowpass filtering and sample rate compressing, 

25. The method of Claim 24, wherein said acts of lowpass filtering and sample 
rate compressing are controlled by an adaptive algorithm. 

26. The method of Claim 20, further comprising the acts of : 

decimating in a first decimator said sampled composite signal before 
providing said sampled composite sigoai to said second iiq>ut of said donodulator, 
wherein said first decimator has a first lowpass filter transfer fimction and a first 
decimation rate; and 

decimating in a second decimator an output of said demodulator, 
wherein said second decimator has a second lowpass filter transfer function and 
a second decimation rate. 

27. The method of Claim 26, finflier comprising the act of controlling said first 
lowpass filter, said second lowpass filter, said first decimation rate, and said second 
decimation rate by an adaptive algorithm. 

28. The method of Claim 26, wherein a product of said first decimation rate 
and said second decimation rate is substantially constant 
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